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CHAPTER 1 
INTRODUCTION 


1 • 1 PIEZOELECTRICITY: 

P lezoelectrici ty is the creation of an electric charge by the 
applied stress. Cady 1113 defined piezoelectrici ty as 'electric 
polarisation produced by mechanical strain in crystals belonging 
to certain classes, the polarisation being proportional to the 
strain and changing sign with it'. Two effects are observed in 
piezoelectricity - the direct and the converse. The direct effect 
IS the creation of an electric charge by the application of stress 
on the piezoelectr 1 C solid. In terms of dielectric displacement 
'D' (charge Q per unit area) and stress 'T' we may write 

D = ^ = dT (d expressed in coulombs/New ton ) 1.1 

The converse effect is manifested when an applied field 'E' 
produces a proportional strain 'S', expansion or contraction 
depending on the polarity. 

S = dE (d expressed in meters/volt) 1.2 

For both effects, the proport lonal i ty constant is the 
piezoelectric constant 'd' which is numerically identical in both 
cases. 


d 


D „ S 
T " E 


1 .3 


Non~cent rosymme t ry is the essential condition for piezoelectricity 
with only one exception. 21 out of the 32 crystal classes are 
non-centrosymmetr 1 C . 20 out of 21 crystal classes show 
piezoelectricity . The single exception is (432) class of 
crystals. A further subgroup can be made, consisting of those ten 
classes which contain a unique polar axis (an electric dipole 
moment) in the unstrained condition. This is the pyroelectric 



subclass, Fig. 1.1 summarises the interrelationship of 
piezaelectr ICS and subgroups on the basis of internal crystal 
symme t ry . 


PIEZOELECTRIC CONSTANTS : 

The piezoelectric constant 'd' which relates the dielectric 
displacement 'D' to the strain 'S' is probably the most 
fundamental for conveying the concept of piezoelectricity. 
Another frequently used piezoelectric constant is 'g' which gives 
the field produced by stress and is known as piezoelectric voltage 
coe f f ic lent , 

_ E _ S r Volt/meter _ Volt /meter 

^ T D [Newton/Sq. Meter Newton 

The 'g' constant is related to the 'd' constant by the 
permi t t i v i ty . 




Ks 


1.5 


'K' IS the relative dielectric constant of the material. 
Additional piezoelectric constants which are only 
occasionally used are 'e' which relates stress 'T' to field 'E', 
and 'h' which relates strain 'S' to field 

T = -eE l.h 

E = - hS 1.7 

Rigorous development of these relationship may be found 
elsewhere C23. Jaffa C43 defined the piezoelectric constants as 
partial derivatives evaluated at constant stress (subscript T), 
constant field (subscript E), constant displacement (Subscript D) 
or constant strain (Subscript S) . These boundary conditions can 
be thought of as 'free', 'short circuit', 'open circuit' and 
'clamped' respectively. The actual definitions are 


f«l - 



N 


d 


1 .8a 




Figure 1,1 : I n ter re 1 t lonshi p of piezoelectrics crys 

subgroups on the basis of internal crys 
symmetry 


ft- 





(-S). 


1.8b 

1.8c 

1 .8d 


The equations of state relating the electric and elastic 
variables can be written in a general form as follows: 


D = dT + £■ E (Direct piezoelectric effect) 1.9 

E 

S = S T + dE (Converse piezoelectric effect) 1.10 

The elastic, dielectric and piezoelectric properties are 
orientation dependent in all piezoelectric media. For the 

symmetry of ceramics (a mm), the general equation is replaced by a 
specific set of equations. 

D. = ^.E. + d_T„ 1 . tla 


effect ) 


■^1^1 '^Ib'^S 


^2 ^ *^ 1^2 ^ 15''^4 

^3 " ^ 3^3 *^31 ^ 33‘’’3 

®i = ®11^T, ®12^T„ ^13^T^ ^ ^31^3 

X jC Of 

S* — S. j E_ S . _E— + S^-.E_ d™ E_- 

2 11 T_ 12 T, 13 T_ 31 3 

2 13 


1.11a 
1.11b 
1.11c 
1 . 12a 

1 . 12b 


= S,_E,_ - . + S_-,Et. + dT—E^j. (Converse effect) 1.12c 

3 13 33 I ^ 33 3 

12 3 


^4 ~ ^44^T^ ’^15^2 

®44^T "^15^1 

o 

S = S T 
6 66 & 


1 . 12d 
1 . 12e 
1 . 12f 


The subscript 3 refers to the poling axis and 1 and 2 refer 
to arbitrarily chosen orthogonal axes in the plane normal to 3. 
Subscript 4, 5 and 6 refer to the shear stress and strain in 
planes normal to the 1, 2, and 3 axes, respectively. Also, by 
conventional notation, the first subscript of the 'd' constant 


gives the electrical direction (field or electric displacement) 



while the second gives the component of mechanical deformation* 

Another constant which gives the measurement of the strength 

of a p lezoelect PIC effect is the electromechanical coupling factor 

k* When an electric field is applied it measures the fraction of 

the electrical energy converted to mechanical energy (or vice 

versa when a crystal or ceramic is stressed)* 

2 

The actual relationship is interms of k . 

^2 _ Electrical energy converted to mechanical energy 
Input electrical energy 


or 

^2 _ Mechanical energy converted to electrical energy 

Input mechanical energy 

As the conversion is always incomplete the value of k is 
always < 1 * Typical values of k are 0*1 for quartz, 0*4 for barium 
titanate ceramic, 0*3 0.7 for Pb <Zr Ti. >0.-, 0<><<i, and as 

high as 0.9 for Rochelle salt at its curie point (24^0. 

Dielectric constant is another important property of these 
materials. The relative dielectric constant, K, or 

ratio between the charge stored on an electrode of this material 
brought to a given voltage and the charge stored on a set of 
identical electrodes separated by vacuum. The relative dielectric 
constant is dimensionless quantity. It depends on the orientation 
of grams of piezoelectric materials and differs widely when 
measured parallel or perpend icul ar to the polar axis Cll. The 
dielectric constant when measured at constant stress is called 
'free' dielectric constant denoted by superscript 'T' and when 
measured under constant strain is known as 'clamped' dielectric 
constant denoted by superscript 'S'. The free and clamped 
dielectric constant may differ greatly for strong p lezoe lec t ric 


materials 


The relation between the two is 



= K^(i - k 1.13 

P 

where is electromechanical coupling coefficient (planar)- 

A convenient method C31 of uniquely specifying all these 
directional properties is to use subscripts which defines 
orientation direction as described in Fig- 1-2- Coefficients that 
relate two variables (e-g. the 'd' constant) are indicated with 
two subscripts, the first of which refers to the electrical 
direction (electric field or dielectric d isp 1 acemen t ) and the 


second to the mechanical direction (stress or strain)* Also when 


a property is being measured holding another quantity constant, a 
superscript symbol is used to indicate the quantity held constant. 
The mostly used piezoelectric coefficients with appropriate 
subscripts and superscripts are given below: 

Stress or strain direction 
Displacement or electric field direction 
(electrodes are perpendicular to the '3' axis) 

E 

S = Compliance at constant electric field (electrodes 

shorted ) . 

= Compliance at constant displacement (open circuit) - 

- Dielectric constant at constant stress (free). 

S 

iS: = Dielectric constant at constant strain (clamped). 

= Relative dielectric constant at constant stress (free). 

S 

K = Relative dielectric constant at constant strain 



(cl amped ) « 

d^^ = Charge or strain coefficient (longitudinal); (Charge 

density/ applied stress, c/N; strain developed / applied 
field, m/V ) . 

d^^ = Charge or strain coefficient (lateral), m/V or c/N. 

2 

S.^ - Mechanical compliance (lateral) m /N. 

w 1 
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= Voltage coefficient (lateral), V-m/N (Open circuit 
f 1 e Id / app 1 1 ed stress). 

K-^ = Relative dielectric constant, longitudinal, no units <K^ 

" >^33 = ^33'''=o' 

= Electromechanical coupling factor (longitudinal ) , no 
units* 

k = Electromechanical coupling factor (lateral), no units* 

= El ect romechan ical coupling factor, planar (radial), no 

units (thin disk only), 

Q = Mechanical quality factor, 

m 

Both static and dynamic methods of measurement are used to 
evaluate the various piezoelectr ic constants, A scheme for such 
measurements can be obtained from Appendix 1. 

The details are provided in C41. Only the resonance method 
has been included here to emphasise its superiority over other 
methods of measurement. 

Simple measurements of resonant and antiresonant frequencies 
are used to evaluate piezoelect ric effects. Elastic bodies show 
numerous resonances, the most pronounced are those where the body 
can just accommodate one half wave length of a standing elastic 
wave. The piezoelectric effect is a convenient way to excite such 
elastic waves, to permit observation of the interaction of the 
mechanical resonance with the electric behaviour. This is best 
seen in the equivalent circuit idealised for a single resonance 
(Fig. 1*3), The values of and C, are such that at the resonant 

frequency ‘fr' the impedances 2nfrL^ and -i/27zfrC^ are opposite in 
sign and equal in magnitude, and the total impedance of this 
branch is given only by the mechanical resistance ' » which is 




Fig. 1.3 (a) Equivalent circuit of a piezoelectric body 
(b) Reactance of piezoelectric resonator. 



comparatively quite small. This mechanical branch is however in 
parallel with the electric capacitance 'Co'. The frequency of 
minimum impedance of * resonance ' , 'fr' is that of a standing wave 
under 'zero field' conditions. For applied frequencies above 'fr' 
the mechanical branch becomes inductive. When the series 
impedance of this branch becomes opposite and equal to the 
impedance of 'Co', the impedance of the crystal reaches a maximum 
since opposite currents flow in the two branches. This conditions 
of parallel resonance, however means constant charge on the 
crystal. The frequency of maximum impedance or ' an t i resonance ' 
'fa', IS therefore that of a standing wave under open-c i rcui t 
conditions. 

1.3 LEAD-ZIRCONATE - TITANATE CPZT> SYSTEM : 

The most important piezoelectric ceramics <PZT) crystallise 
in the Perovskite structure. This structure may be descr ibedl as a 


simple cubic unit cell with a 

large cation 

(A) 

CPb'*'^, 

Zr 3 on 

the 

corners, a smaller cation <b) 

in 

the 

body 

center , 

and 


oxygen (O) in the centres of the faces CFig. 1.43. The structure 

is a network of corner linked oxygen octahedra with the smaller 

cation filling the octahedral holes and the larger cation filling 

the dodecahedral holes C43. A wide variety of cations can be 

substituted in the perovskite structure C51. The relationship 

t = (R 4- R )/ yz (R^+R^) 1.14 

a Q BO 

describes the ideal perovskite structure where t = 1.0 and R * , R^ 

JH O 

and Rq indicate the ionic radius of the large cation, small cation 
and anion respec t i ve ly . In practice those structures whose 
tolerance factor t = 0.95 to 1.0 are cubic; those with lower t 
values are slightly distorted but nonf erroe lect r ic . 




Polansabil ity of the ions is also an influential parameters. 
While considering this 'hardsphere' model of ion in contact, it is 
also necessary to consider that the bonding is not necessarily 
completely ionic. 

The discovery of very strong and stable p lezoelectr ic effect 
in PIT solid solution L61 in 1954 gave a jolt to the researchers 
in this field. Higher curie temperature of pZT <360^0 over 
BaTiO-y <120'^C) gave stable piezoelectr ic properties over a wide 
range of temperature. 

PZT IS a solid solution of lead~zi rconate (rhombohedral ) and 

lead-t itanate (tetragonal) with varying degree of composition. 

The generalised chemical formula is PbCZr Ti^ ) where 0 <k<1- 

K l-x 3 

Esturo Sawaguchi 171 first developed the phase diagram of the 

solid solution of PbZrO_ - PbTiO_ system. He established that 

•+ "+*4 

substitution of Zr for Ti in PbTiO^ reduce the tetragonal 

distortion CFig. 1.53 and ultimately causes the appearances of 
another ferroelectric phase of rhombohedral <R.3m) symmetry. The 
boundary between tetragonal and rhombohedral phases is nearly 
independent of temperature (morphot rop ic ) . The current version of 
phase diagram CFig. 1.63 is determined with zirconia containing 
low amount of hafnia. 

P. Ari-Cur and L. Benguigui C83 showed by means of X-ray 
diffraction the co-existance region to be a wide range of 
composition. Initially, it was thought that tetragonal - 

rhombohedral phase boundary (MPB) corresponds to a specific 
composition. 

The solid state reactions leading to the formation of PZT 
solid solution have been studied by many investigators but with 
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Fig. 1.5 Unit cell distortion at room temperature tor 
the PbTi03 -PbZr03 system. 




Fig. 1.6 PbTi03 ~ PbZr03 sub-solidus phase diagra 




differing conclusions. According to Matsuo and 


Sasaki C93 when 
PbO, T 1 O 2 and ZrO^ in the molar ratio 2:1:1 is heated; PbO and 

TiQ_ react to form PbTiO_ almost completely at 650^C. The 

Z. o 

raaction PbTiO-y + PbO + ZrO^ > Pb (Zr Ti, ) 0-. starts near 

32 K l~x 3 

650'^C« In this reactions soon as ZnO^ starts to react with PbO, 

PbTiO^ reacts to form Pb<Zr Ti. )0^ solid solution, 

3 X 1-K 3 

The PZT ceramics offer good piezoelectric properties over a 
wide range of compos i t ions . Jaffe et al C63 reported the 

excellent performance of PZT ceramics with compositions near the 
MPB. They observed high value of dielectric constant and coupling 
coefficient as well as relatively small variation of these with 
changing temperature. Composition which lies in the MPB region is 
of particular commercial interest, 

1,4 EFFECTS OF IMPURITY DOPING IN PZT CERAMICS : 

Lead zirconate titanate, Pb<Zr, 11)0^, ceramics are excellent 
piezoelect r 1 C materials. Their properties can be improved by 
doping with a small quantity of additive oxides. These dopant 

atoms either replace partially the constituents atoms or they make 
pseudoternary solid solution compounds. Properties can be further 
tailored with multiple oxide dopings. 

B. Jaffe C43 suggested that most of impurity atoms enter into 
the lattice points to a certain extent and thus the ceramic 
piezoelectric properties change. Depending on the substitution of 
cations in ABO^ perovskite type structure of PZT, the dopants or 
modifiers have been classified as, 

1 . Insolvent cation substitution type 

2-i- 2+ 2"^ 

Example : A site (Pb ) substitute by Sr , Ba 


11 . Donor type 



M 4 . 

Mo etc 


Example A site Bi^^,Th^^ 

B site ; Sb^-^, Ta^-", , 

111 . Acceptor type 

A site ; k"*^, Na"*^ 

B site : Mg^-^, Sc^'" , Fe^^, Al^"-, Cr^'" etc 
Ion valences of the acceptor atoms are lower than those of 
the constituent atoms. Thus oxygen vacancies are introduced for 
valence compensation. Ceramics doped with acceptor atoms have 
fol lowing propert les: 

1, Relatively low dielectric constant, 

2, Low dielectric loss, 

3, Moderately lowered electrical resistivity, 

4, High mechanical quality factor, 

5, High coercive field. 

6, More difficult poling and depoling, 

7, Relatively dark color, 

8, Relative insensitivity to darkening by Alight, 

Donors produce variances in Pb~sites, because their ion 
valences are higher than those of the constituent atoms. Ceramics 
doped with donors have a different set of properties that are 

1, Increased dielectric constant, 

2, High dielectric loss, 

3* X3!lj€jiied elastic compliance, 

4, Low mechanical Q, 

5, High piezoelectric coupling factor, 

6, Low coercive fields relatively square hysteresis loops, 

7, Greatly increased electrical volume resistivity. 

Anomalously small aging effects. 
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9. Easy nan elastic mechanical deformation* 

10* Yellow color* 

11. Translucency * 

12. Easy phototropic darkening. 

The acceptor substituents thus have opposite effects fron 
those of the donor substituents. 

It IS difficult to assign valencies to Cr and U additives C43 
in Pb(Ti, Zr> O ceramics. The properties suggest that these 
atoms are present in the lattice in more than one valence state* 
The properties of the ceramics are partially those of the donor 
doped and partially those of the acceptor doped ceramics* The 
following charac t er ist ICS are noted for Uranium doped PZT 
ceramics , 

1. Low aging (sometimes negative) 

2* High resistivity at elevated temperature* 

3* Slightly more conductive at room temperature than unmodified 
compos 1 1 ion * 

4* Mechanical Q very high for rhombohedral compositions. 

5. Dissipation factor (IKHz) high. 

6* Coupling factor normal. 

and for chromium doped the properties are, 

1* Low aging. 

2. Resistivity slightly lowered. 

3* Mechanical Q slightly high. 

4. Dissipation factor (IKHz) high* 

5. Reduced coupling factor. 

6. Rslitivsly lo« dependence of dn tempereture for 


selected compos i t ions* 



Many of the properties of donor doped PIT materials ar 
explained on the basis of A site vacancies facilitating domai 
boundary motion. Domain walls can be displaced with very smal 
electric field or mechanical stresses. This increases bot 
elastic and dielectric compliance and loss factor^ The easy wal 
motion greatly facilitates initial stress relief. This explain- 
the low aging rates, since those domains under stress switch t< 
relieve most of the stress immediately upon removal of the field. 
The case of switching is related to relatively low coercive fielc 
and squareness of hysteresis loops. 

In the case of acceptor doped PIT materials vacancies in the 
Oxygen position will tend to shrink the cell size. This shrinkage 
and distortion of the cells by Oxygen vacancies is responsible for 
increasing the mechanical quality factor, high coercive force and 
for lowering dielectric constant. 

Thus domain wall motions control many important properties. 
Poling process takes place owing to the domain wall motions. 
Mechanical quality factor depends on the domain wall hysteresis 
loss. Aging effect is also due to the gradual domain wall 
motions. Internal bias field which is controlled by impurity 
doping, has an influence on the domain wall motions. 

An asymmetric DE hysteresis loop in poled BaTiO^ ceramics is 
observed by investigators and they proposed the existence of an 
internal bias field, presumably arising during the aging process. 
Subsequently the same phenomenon has been found in Pb (Zr,Ti)0-^ 
ceramics by Takahashi C103. Properties like mechanical quality 
factor and the aging change in frequency constant and also 
dielectric constant are affected by the internal bias field C 1 1 1 



controlled with various single impurity additives in Pb(Zr,Ti)0^ 
ceramics* 

One of the most important properties of piezoelectric 

resonators is its mechanical quality factor* Mechanical quality 

factor depends on the various impurities present in lead zirconate 

ceramics* Qm is defined as the reciprocal of internal friction 

C123. Acceptor Oxides such as Fe^O-^, Cr^O^ and ^^^^2 

Pb(Zr, Ti)0-^ ceramics in a small amount, were found to increase 

the values of Qm, However, donor oxides such as La_0 , Nb_Q-, 

2 3 2 b 

ThO^, WO^ , and Sb-0- decrease the values of Qm* 

Decrease of Qm can be explained by the following proposed 
mechanism. Impurities, which create vacancies in Pb position 
reduces the local stress in the domains which can therefore 
undergo domain switching easily* As a result the domain walls 
become remarkably mobile and produce only micro-hys te res is * It 

also suggests that Qm decreases with increase of internal friction 
and hysteresis. 

In the other case, that of acceptor doping, increase of Qm 

can be explained by the presence of large amount of space charge, 

most probably due to the oxygen vacancies, which suppresses the 

domain wall motion and decreases the internal friction. 

The amount of space charge can be estimated ClOl by an 

observable parameter which is defined to be the reduced amount of 

space charge (P - P )/P , where P = the value of saturated 
sis s 

polarisation and = the value of polarization determined by the 
hysteresis curve obtained just after the application of an a.c. 
field* The values of <P - P . ) /P . vary with the kind and amount 
of impurities* Poling efficiency was found to decrease with the 



increase of <P -P ) /P and the lowering of poling temperature. 

Si s 

Qm can be determined by means of the piezoelectric resonance 
With a radial fundamental mode* In accordance with appropriate 
IRE standards, the values of Qm can be determined by the 
measurements of resonance frequency fr, an t i resonance frequency 
fa, resistance at resonance R and capacitance C* Thus values of 


Qm IS obtained by means of the relation given below, 



2nf r (fa^-fr^)RC 


where the value of C is determined at IKHz* 


1*15 


Specific resistivity (p) has been measured applying a dc 
field of lOoV/mm to ceramics* Sadayuki Takahashi and H* Takahashx 
have shown C121 the variation of mechanical quality factor Qm of 
PbCZr^ csramics containing various kinds and amounts 
of impurities with a reduced amount of space charge • 

They found that the values of Qm 
increases linearly with (P -P ) /P at constant specific 

SIS 

resistivity* These are mainly applicable to the impurities of 
acceptor type* 

Electromechanical properties of piezoelectric materials can 
be improved by doping with impurities in order to meet the 
requirements for particular applications. It is said that some of 
the doping effects on the electromechanical properties are due to 
an internal bias field controlled by impurity doping* S. 
Takahashi has studied the impurity doping effects on Pb (ZrTi>0-^ 
ceramics from the internal bias field view point* 

Ferroelectr 1 C DE hysteresis loop is symmetric with respect to 
the origin for usual ferroelectric as shown in Fig 1*7 (a). If a 
d*c. bias field is super imposed on a measuring a*c* field, the 



loop IS displaced asymme t r ical ly with respect to the origin by E 

as shown in Fig l,7(b) below. The loop becomes like that of Fig 

1.7(b) for some f erroe lec t r ic ceramics, although no d.c. bia 

field has been superimposed. It appears as if the d.c. bias fiel 

exists naturally inside of such ceramics. This field is definei 

as the internal bias field E . Thus the E value can bj 

1 1 

determined from the hysteresis loop. 

l.S IMPURITY DOPING EFFECTS ON INTERNAL BIAS FIELD 

According to S. Takahashi Cli3 internal bias is strengthenec 
by doping with the following impurity atoms, Fe, Ni , Co, Ga, In, 
Cr, Rh , Ir or U, but weakened by La, Mb, Ta, Sb , Bi, Th or W. 
Among these atoms, Fe, Ni,Co, Ga, and In are known as acceptors 
and those which weakened the internal bias field are known as 
donors. Rh and Ir atoms behave like Cr and U atoms present in the 
ceramics. S. Takahashi C113 selected the impurity atom K, Fe for 
acceptor group, La, Nta for donor group and Cr for remainderCxe 
partially acceptor and partially donor) and examined the doping 
effects on the internal bias field. 

The internal bias field (Ei) value increases with increasing 
doping concentration (x) for K or Fe while it decreases slightly 
for La and Nb . For Cr doped ceramics E^ increases with x when x 
IS less than 0.5 at %, but decreases when x is more than 0.5 at %. 
This decrease may be due to the internal bias field relaxation 
during the course of switching. In the case of multiple impurity 
doping it was found that doping with same kind of dopant gave rise 
to effects similar to those in single impurity doping. Also 
equimolar mixed doping of acceptor and donor seems to exert no 


influence on the Ei value. 




1.6 INTERNAL BIAS FIELD ORIGIN 


Lambeck' and Jonker [13,141 found that, in reduced Hn-doped 
+2 

BaT jlO^ crystals, Mn xons substituted for Ti-sites and Oxygen 
vacancies form complex with dipole character. In this case 
dipoles gradually orient themselves parallel to the spontaneous 
polarisation Ps and thus they stabilizes the polarisation. For 
PZT ceramics, acceptor doped atoms create, oxygen vacancies by 
valence compensation. Thus as in the case of BaTiO^ doped with 

O’ 

+2 

Mn , here complexes of acceptor ions and induced vacancies 
stabilize the spontaneous polarisation and are the origin of 
internal bias field. 

But according to S. Takahashi Cill, Oxygen vacancies 
introduced by acceptor doping move slowly to the domain walls by 
diffusion, they pm the walls and stabilize the domains. This 
also leads to asymmetry in the switching current. At present it 
IS certain that only acceptor doping is responsible for the 
generation of the internal bias field. Experiment shows that 
donors are ineffective in generating an internal bias field and 
are also ineffective in stabilizing the polarization. This is 
because they produce Pb~sites vacancies, which facilitates motion 
of the domain walls under electrical or mechanical drive. 
Internal bias field is little changed in the case of doping with 
equal amounts of acceptor and donor Ci51 from that in undoped 
ceramics . 

Impurity doping also effects the poling process. This was 
examined from the internal bias field view point. Thus from the 
DE hysteresis loop domain wall motion under a strong electric 
field can be found out. DE hysteresis loops for modified PZT 



ceramics doped with a single impurity like 5 atV* K, 0.8 att4 Fe, 
2.4 at%La, 2.4 at% ^4b and 0.53 at% Cr and also loop far undoped 
PZT are shown in Fig. 1.8. The loops were observed at room 
temperature by applying a 50 HZ a.c. field with 4 KV/mm peak 
value. The undoped and donors doped PZT show saturated normal 
loops just after applying the field, but the acceptors doped PZT 
show unsaturated constricted loops for a long time. For Cr-doped 
PZT, a constricted loop is obtained just after applying the field. 
Then, the loop changes into a saturated normal one in a short 
time. The unsaturated constricted loop may be due to the 
existence of the internal bias field. 

The change in the loop during the application of the field is 
due to the internal bias field relaxation. Relaxation time 
constant for the internal bias field is smaller for Cr-doped 
ceramics than for acceptor doped ceramics. Thus from this it can 
be expected that because of the strengthened internal bias field, 
domain wall movement is difficult in the acceptor doped PZT. 
Domain walls move easily in undoped and donor doped PZT even by a 
d.c. field application. As poling is carried out accompanying the 
domain wall motions so internal bias field effects on poling 
process. 

Also mixed impurity doping effects on mechanical quality 
factor Qm and aging change in frequency constant Nr. The 
experimental results £1163 of doping effects on the properties of 


ceramics are summarized in Table 1.1s 




ceramics doped with single impurity 


Table 1.1 s Doping effects on ti*e elect romecbani cal properties 


of ceramics 


Dopant 

E 

1 

m 

iiNr/Nr 

p 

T 

^33 

K 

r 

a 

I 

I 

I 

D 

ID 

ID 

d 

D 

D 

D 

I 

ID 

ID 

a+a 

I 

I 

I 

D 

ID 

ID 

d*^d 

D 

D 

D 

I 

ID 

ID 

a+d (a>d) 

I 

I 

I 

D 

ID 

ID 

a-i-d (a==d) 

UC 

UC 

UC 

UC 

ID 

ID 

a*+-d ( a<d ) 

D 

D 

D 

I 

ID 

ID 

PFN 

UC 

UC 

UC 

UC 

ID 

ID 

PFN+a 

I 

I 

I 

D 

ID 

ID 

PFN+d 

D 

D 

D 

I 

ID 

ID 

where Base 

Dopant : a 

material : Pb 

: acceptor <K 

(2r 

or 

0.52'^^0.4e^*^3 
Fe ) ^ d : donar 

< La or 

Nb> 


PFN : Pb (F.i,2Nb,^2)03 

Symbol for doping effect 

: I 

: increase; 

D ; 

decrease 

*, UC. : 

unchange ; 

It IS 

ID : increase i 

obvious that 

□r decrease. 

Qm and the aging change are 

closely 


correlated with the internal bias field E . As for the other 

1 

quantities there seems to be no correspondence with E^. 

In Ferroelectric ceramics, the mechanical loss Qm^ is mainly 
due to dielectric loss, so Qm is linearly related to tan <5 through 
electromechanical coupling- On the other hand, dielectric tan <5 
IS the sum of macro-hysteresis loss, micro-hysteresis loss and 
intrinsic dielectric loss. The internal bias field reduces the 
hysteresis loss, through prohibiting domain wall motions. 





1.7 LOW AMPLITUDE APPLICATIONS 


Before discussing the principles of piezoelectric materials 

in electric wave filters certain points in terminology need 

clarifying. According to J.li. Herbert L171 'band width' of a 

filter IS the range of frequencies it will allow to pass with 

minimal attenuation. For a resonating body it is the difference 

between the frequency giving maximum reactance ( an t i resonance > and 

that giving minimum reactance (resonance) and is proportional to 
2 

k . The sharpness of the impedance/f requency peak in the 
neighborhood of resonance is controlled by the losses in the 
material. The lasses are due to both mechanical and electrical 
causes but mechanical loss is taken as the dominant factor near 
resonance and is usually expressed by its inverse, Qm <Q=l/tan <5) . 
Qm can be determined from the width, Af , of the 
impedance/f requency peak at an impedance level which is equal to 
the maximum impedance divided by root two, i.e. Qm = fr/Af where 
fr IS the resonant frequency. Thus k and Qm are the important 
parameters determining the impedance/ frequency charac ter ist ics of 
resonating systems and the term band width may apply, in different 
senses to both of them. In what follows ‘band width' is the 
frequency interval associated with k. 

Solid bodies can be excited into a wide range of vibrational 
modes by suitable excitation 11173. If the excitation is provided 
by a piezoelectric transducer energy will be transferred to the 
body most effectively at its resonant frequency. A second 
transducer may be used to convert the mechanical vibration back 
into an electrical signal and the transfer of the signal from one 
transducer to the other will only occur efficiently at the 



resonant frequencies of the intermediate body* Man> 

electro-mechanical filters have been based on this principle using 
materials with low temperature coefficients of their elastic 
constant (e.g« the Invar steels) for the vibrating bodies and so 
shaping them that only single resonances occur in the waveband of 
interest. Above 1 MHz the resonant condition requires dimensions 
below 1 mm for most materials and, since this dimension needs to 
be controlled with the same precision as that required for the 
centre frequency of the filter, a considerable fabrication problem 
arises, though this can be alleviated to some extent by using the 
natural resonant frequency as a means of monitoring the relevant 
dimension . 

Quartz was the first material in which the mechanical and 

electrical functions were combined Ci73 : the natural crystals 

provided a stable, reproducible raw material with, for specific 

cuts, near to zero temperature coefficients- Advantage was taken 

4 6 

of its very low mechanical losses (high mechanical Q = 10 to 10 ) 
by making filters with sharply defined frequency boundaries, while 
its low piezoelectric coupling coefficient was compensated by 
combining crystals with inductor — capacitor networks when wide 
band width devices were required- 

The far higher piezoelectric coefficients of f erroe lectr ic 

materials led to the expectation that filters could be devised 

with a single body giving the entire filter character ist ics since 

2 

band width is proportional to k <k = e 1 ec t romechan ical coupling 
coefficient). This has indeed been the case for some wide 
tolerance applications but the lower mechanical Qs of ceramic 
bodies and the higher temperature coefficients of 


elastic 



constants of both single crystals and ceramics has led to 
limitations in high stability, close tolerance devices. 

Now the following section £173 describes typical applications 
using plates of p lezoe lectr ic ceramics (0.1 - 1 MHz): 

Discs which are fully silvered on both faces and poled 
normally to their electrodes have radial modes for their lowest 
frequency resonances. The frequency factor, Nr, for the lowest 
radial mode of PZT is about 1 KHzm so that discs of 10 to 2 mm 
radius can cover the frequency range 100 to 500 KHz. The 
fundamental resonance for the thickness mode occurs at a frequency 
higher by a factor of approximately r/t, where r is the radius and 
t the thickness of a disc. 

If the reactance of a disc is plotted against frequency the 
charac t e r ist ic around resonance is as shown in Fig. 1.3. 

Now equivalent circuit for the piezoe 1 ect r ic elements is 
shown in Fig. 1.3. 

Here R^, L^and represent the mechanical losses, inertia 
and elasticity and Co represents the electrical capacitance of the 
element . 



According to Tetsuro Tanaka £183 among piezoelectric filters. 


there are quartz filters, 


which have the longest history and 



comparat 1 ve ly new filters such as mechanical filters, ceramic 
filters and SAW (Surface Acoustic Wave) filters. Although quartz 
filters are of the highest quality because they have a high Q and 
excellent frequency temperature character ist xcs , they are used 
only for narrow band filters because they have a small 
electromechanical coupling coefficient and they are expensive, 
LaTiO-^ and LiNbO^ instead, can be used for wide band filters 
though they are also expensive. Quartz filters are mainly used 
for commun icat ion in a high frequency range. 

M. Onoe C193 suggested that mechanical filters utilize a high 
Q and the excellent temperature charac ter ist ics of a constant 
elastic modulus alloy in combination with piezoelectric ceramics. 
They are of high quality but are not suitable structurally for 
very high frequency filters. They are rather suitable for low 
frequency filters, and they can offer even several hundred Hz, for 
example, by using the structure of a tuning fork. 

According to Tanaka C183 ceramic filters are made solely of 
p 1 ezoe lec t r 1C ceramics, whose Q value and temperature 
character ist ics directly determine their quality. They are 
suitable for filters used in consumer electronics and are produced 
in greater number than all other filters because they have a 
considerably wide frequency range and because they are compact, 
strong and inexpensive, 

SAW filters are the newest piezoelectric filters and utilize 
elastic waves which propagate on the surface of piezoelectric 
materials. Their frequency is determined by the dimensions of 
interdigital electrodes, whereas that of the other piezoelectric 
filters are determined by the shape and dimensions of the 



material. Their quality is, therefore, determined both by the 
charac t er ist ICS of the surface of p i ezoel ec t r ic materials and by 
the dimensions of electrodes. SAW filters are suitable 

principally for high frequency range and higher than GHz frequency 
ones are possible. They have a wide applicability ranging from a 
delay line to information processers. 

1.9 PRINCIPLES OF PIEZOELECTRIC CERAMICS FOR MECHANICAL FILTERS 

Piezoelectric ceramics are prepared for fabricating the 
electromechanical transducers used in the mechanical frequency 
filters that find application in long-habul communication systems. 
Also Mechanical frequency filters can offer additional advantages 
such as superior stability, smaller bulk and lower fabrication 
cost. Mechanical filters of uniformly high quality can be applied 
to the automatic mechanics, usually computer controlled, also for 
the fabrication of the filters and component elements, as well as 
for frequency tuning etc. 

First filter developed was 50 KHz channel filter 
specifications C203 of the materials for that is listed in the 
Table 1.2. 

Table 1.2 s Specifications of Piezo-ceramics for SO KHz — filter 

Permittivity s/e % 900±40 

Coupling coefficient = 0.25 ± 0.02 

Figure of merit k^, = 56 

£4 X 

o 


Resonance Q 

Temperature Coeff, of 
Permittivity 

Temperature Caeff^ of 
Frequency Const an t 


Q >> 600 


Tc^ 

TcF 


^ (60'^C) - fr(20'^C) ^ 

si20°C) 40K 
fs<6Q°C) - fs<20'^C> 

fs<20°C) .40K 


3000 ppm./K 
= 80±50 ppm/K 


Aging of frequency AF<0«iy4 

Constantper time decade 

The resonance Q of the ceramics for the new signal filters 



has to be at least three times higher. Also the temperature 

coefficient of the resonance Q, TcQ, must be very small. So for 
fabricating filter material we have to dope pure PIT with dopants 
like Fe, Cr, Mn or other acceptors which result high mechanical 
resonance Q. Also the high stability with respect to time of the 
frequency constant of ceramics doped with Cr or Hn makes them of 
interest for application in frequency filters. For resonators and 
filters temperature coefficient of frequency constant (TcF) values 
ranging between 0 and 80 ppm/K are required. 

Principal properties of Mn doped Pb (2r Ti^ )0_ ceramics for 

K l~x 3 

filters, as well as a portion of the phase diagram of these 
ceramics are discussed by H. Thomann and W. Wersing C201. As with 
all lead zirconate titanate ceramics, the dielectric, 
p 1 8 208 1 ectr 1 C and elastic properties and their temperature 
coefficients are closely related to the various f er roe lect ric 
phase transitions. The morphotropic tetragonal /rhombohedral phase 
boundary is of special interest because the permittivity and 
electromechanical coupling coefficient of the ceramics here reach 
their peak values. This is because spontaneous polarisation in 
the rhombohedral phase with its 8 possible directions can take on 
a more favorable orientation with respect to the applied electric 
field than is passible in the tetragonal phase with its 6 
directions and also because the spontaneous strain in rhombohedral 
materials is less than in tetragonal, the rhombohedral phase is 
both electrically and elastically privileged during the 
polarization of such ceramics. 

Ceramics doped with Cr, Mn and V give stabilizing effect. 
But for these, ceramics near the morphotropic phase boundary would 



be generally unsuitable for filter application* Phas) 

stabilization is caaaally related to the stabilization of th^ 
domain configuration* The s t ab i 1 i zat ion of the domair 

configuration in ceramics is directly indicated by their 
hysteresis loops* Thus non polarised samples exhibit constricted 
and polarized samples exhibit displaced hysteresis loops* This 
effect IS usually described by introducing an internal stabilizing 

field E which is discussed earlier. 

1 

All symptoms of the stabilizing effect that is linearization 
of the TcF, reduction of aging rate, appearance of an internal 
field appear to have the same physical origin- 

^4ume^aus models were discussed in the literature for the 
purpose of explaining the stabilization are classifiable under 
three basic concepts: 

1. Bulk Effect 1211 

Lattice defects of a dipole character <8.g. Hn" Ti - 
assume orientation in the Lorentz field of spontaneous 

polarisation and freeze the direction of pol ar isat ion . 
a* Domain Wall Effect C21, 222 

Lattice defects diffuse to the domain walls, there by 
reducing the domain wall energy and freezing the walls* 

3. Grain Boundary Effect tlO, 211 

Charged impurity atoms or vacancies accumulate at 
nonsaturated domain configurations of neighboring grains and 
freeze the direction of polarisat ion* 

It IS obvious that oxygen vacancies appear to be of critical 
significance in all these models and stabilization can not be 
attributed to a single effect in the case of all materials* 



After treating the basic phenomena of significance in 
preparing filter ceramics based on lead zirconate titanate, we now 
discuss how the properties of filter ceramics can be improved by 
the combination of additions. 


For filter ceramics the ternary system like PbTiO -PbZrO- ” 
Pb <Mg ^ ^ 0-^ 9 and many other solid solutions like 

1 /3S'=2y3'°3 

Pb(2r,Ti)0^ are often additionally doped with a stabilizer such as 
Mn . The ternary PbTiO^ - PbZrO^-Pb (Mn , ._Nb_ >0^, system, reported 
for the first time by Wang C231 has composition with favorable 
planar electromechanical coupling factor kp and very high Qm and 
IS recommended for filter ceramics. Also the temperature 
coefficient of resonance frequency <TcF) of the piezoelectric 
PbC(Zr Ti ) (lin Nb^ ^ 30^ ceramics is investigated 
{1233 by ways of Zr/Ti ratio, additive, sintering and poling 
temperatures . The morphotropic phase boundary (MPB) locates at x 
value near 0.51, where the planar electromechanical coupling 
factor (kp) reaches a maximum with TcF value of approximately O 
ppm/*^C* The ceramic composition e’^hibiting TcF values of O ppm/'^C 
IS changed with processing conditions and be explained by the 
shift of the MPB compositions or the phase change. 

Ceramics with optimally high permittivity and 
electromechanical coupling on the one hand, and high resonance Q 
and low aging rates on the other, are obtained by choosing an 


Mn/Nb ratio of ^ 0.2, 


Ceramics with particularly favourable 


reproducibil i ty with respect to the TcF are realized with 


zirconium mole fractions between 0^43 and 0.44. 


Mechanical losses of filter ceramics should be as low as 



possible and the mechanical resonance Q should be as high as 
possible in order to minimize the insertion lossm 

To keep the insertion loss of narrow band signal filters 
sufficiently independent of the temperature, the temperature 
coefficient of the ceramics Q should be as low as possible 
(jTCQ/Q| <5 ppm/K). Mechanical and electrical losses in soft 
piezoceramics mainly arise due to vibration of domain walls and 
their interaction with lattice point defects. Both Serthsen, 
Hardtl and Schmidb and Arlt and Dedenchs have developed models 
that describe the losses as a result of the damped motion of 90"^ 
domain wails. These models lead to the following relationship 

C243 between electrical < tan <5) and mechanical losses (1/Q) 

2 T 

S . 

l/Q = tan a m 1.17 

S = Spontaneous strain 


aa 

= Spontaneous polarisation 
= Permittivity 

OnJ 

F 

= Compliance 

Mechanical losses in piezoelectric materials is related to 
electrical conductivity (.€>) by the equation, 

tan <5 = 1.18 

The electrical conductivity of Mn doped ceramics 

_ 1 Q _9 

(10 -10 l/f2m) IS far higher than that of ceramics doped with 

donors or complex oxides. Mn sometimes shows anomalously large 

dielectric losses due to inhomogeneous high conductivities 
-5 —4 

(10 -10 i/Osn) at grain boundaries. The resonance Q measured on 


such ceramics at about 100 KHz is only around 400. 



CHAPTER 2 

STATEMENT OF THE PROBLEM 

One of the most important properties of piezoelectr ic 
resonators is its mechanical quality factor (Qm>m Mechanical 
quality factor depends on the various impurities present in Leac 
Zirconate Titanate Ceramics* Based on literature survey we have 
chosen four dopants namely , MnO- , Al-0«, 

Using design of experiments methodology an attempt has been 
made in present work to find out the optimum values of the 
different parameters (sintering time, sintering temperature, 
dopants concen t rat ion > needed for a tuner material. 



CHAPTER 3 

EXPERIMENTAL PROCEDURE 


3.1 RAW MATERIALS: 

The raw materials used for the preparation of PZT sample? 
were as follows. 


Table 3.1 : Details of CbewdLcals Used 



Chemicals 

Pur 1 ty 

Manuf ac turer 

1 . 

Lead Oxide (Yellow) (PbO) 

99.9*/. 

Aldrich chemical 
Company, Inc* 

2. 

Zipconia <ZrO^) 

99 . 9*/, 

Magnesium Electrt 
Ltd (MED Grade » 

3. 

Titania (TiO^) 

A. 

99.9*/. 

Aldrich chemical 
Company, Inc. 

4. 

Niobium PentoKide (Nb^Q- ) 

2 5 

99 . 9*/, 

Aldrich chemical 
Company, Inc* 

5. 

Manganese dioxide (MnO^) 

99 . 9*/, 

Aldrich chemical 
Company, Inc. 

6 . 

Bismuth Oxide (Bi^Q.^) 

99 . 9*/. 

Aldrich chemical 
Company, Inc* 

7. 

Alumina (Al^O..^) 

99 . 9*/. 

Aldrich chemical 
Company, Inc* 


3.3 PREPARATION OF ATMOSPHERE POWDER s 

The composition of the atmosphere powder was selected as lead 
zirconate (PZ) (1 mole PbO ; 1 mole ZrO^) + 10 wt*4 PbO. 

250 gms of batch material for the atmosphere powder was 
prepared by weighing PbO and ZrO.^ accurately. The batch material 
was put in a plastic jar. It was wet ball milled with Al^O^ balls 
for 6 hours. Triple distilled water was used as liquid media. It 
was dried in a breaker with constant stirring at 80‘^C ± 5°C, The 
dried powder was poured in a recrystallised alumina crucible and 
covered with a similar crucible. The whole assembly was 





transferred to a pit furnace with silicon carbide heating 
elements. It was calcined at 860^C for 1 hour. After ta^'ing out 
from the furnace the atmosphere powder which formed a compact was 
crushed to powder by agate mortar and pestle* A fast XRD scan of 
this powder confirmed the presence of lead zirconate, Lead Oxide, 
and zirconia phases only. The prepared powder was thus ready to 
be used as the atmosphere powder during the calcination and 
sintering of the PIT samples. 

3.3 PREPARATION OF PZT COMPOSITIONS s 

In the present investigation an attempt was made to obtain 

the optimum values of different variables in order to get maximum 

quality factor Qm and dielectric constant (before poling) of the 
sample. Effect of different variables like sintering time, 
sintering temperature, concentration (atom fraction) of dopants on 
dieiiectric constant and mechanical quality factor were studied. 

3.4 RESPONSE SURFACE METHODOLOGY : 

Response Surface Methodology (RSM) was used for this 

optimization study. This technique has been discussed elaborately 
by Kitttrell and Erjavee C253, Box et al. C261 , Davies C273, Khun 
and Cornell C283 . 

Let Y be the response of a chemical process dependent on the 
levels of K factors X^, which can be precisely measured 

and controlled. The model for the u^^ combination of factor 

levels is given by, 

Y = fx, , X^ , . , Xj. 1 + ^ 

u ^ lu’ 2u’ ’ Kuj u 

u “ 1, 2,...., N 

where, H is the number of experiments 



IS the functional rel at lonship ^ and 


€; IS the error involved. 

A geometrical portrait of the response function in the factor 
space IS called a response surface. The experimental region R, is 
a bounded subspace of the whole factor space. The experimental 
region is bounded because of the practical limitations. The 

response surface methodology locates a point <xf, X^) 

Within the experimental region R, at which V is an extremum 
[29.30.311. The traditional one factor at a time method^ where 
the experiments are conducted holding all the other factors 
constant often fails to locate the true optimum. 

The response surface methodology cuts down the experimental 
effort by making use of experimental designs which permit the 
experimenter to assess the strength of the interactions between 
the factors while varying them simul taneously [301. 

RSM consists at several steps. 

1. The first step is to design a set of experiments and conduct 
them to get reliable estimates of the parameters. When the 
number of factors is small, factorial design can be used. On 
the other hand, fractional factorial designs are used when 
the number of factors are large. The details of factorial 
and fractional factorial designs are discussed by Box et al. 
C263 and Rao and Iyengar C321. 

2. A suitable mathematical model is proposed to fit for the 
experimental data and then test for model adequacy through 
Lack“af““fit F-tests [331. 

3. To find out the optimum conditions of the independent 



variables which will produce the maximum (or minimum) value 
of the response. At a point which is remote from the 
optimum, there is a little curvature in the true response 
surface and first order models will be satisfactory to 
describe the response surface. In the vicinity of the 
optimum, however, higher order models will be required due to 
the presence of curvature in the response surface. 

RSH is sequential in nature. Some experiments are carried 
out, valuable information is gathered and the next stage is 
designed for getting better values of the response. Generally the 
method of steepest ascent (or descent) is applied for moving 
sequentially along the direction of maximum (or minimum) increase 
(or decrease) in response. 

3. 3. a VARIABLES IDENTIFICATION AND THEIR LEVELS : 

Important variables which influence the dielectric constant 
and mechanical quality factor are sintering time, sintering 
temperature, atom fraction of various dopants. The dopants 
selected in the study are Nb-0-, HnO , A1 0^ and Bi 0^. After 
knowing the important variables influencing the process, a base 
level had to be chosen within the experimental region. The base 
levels of the factors were chosen on the basis of an a priori 
knowledge. Further the variation interval for each factor had to 
be chosen. The variation interval of a factor when added or 
substrated from the base levels gives the upper or the lower level 
of the factor respecti vely . 

To simplify the recording of the conditions of an experiment, 


and the processing of the experimental data the values of the 



factors were coded as follows 


X - X 

X = -i__ L2. 1 = 1, 2, , K 

1 

where X is the coded value of the i-th factor 

i 

X IS the natural value of the i~th factor, 

X IS the natural value of the base level of the i-th 

IQ 

factor, and is the variation interval of the i-th factor^ 

Conventionally, the coded value of upper level corresponds to 
•4-1, the lower level to -1 and the base level to O. The upper, 

lower and base levels of the variables in this first stage of 
design are given in Table 3,2. 

The following dependent variables (responses) were considered 
to be important for optimization. 

a. dielectric constant (before poling) 

b. mechanical quality factor. 

Table 3.2 2 First Order Response Surface Strategy CFirst Move> 

LEVEL OF FACTORS 


Factor 

Code 

Base 

level 

Lowe r 
Level 
(-) 

Higher 

Level 

(+) 

Sintering time (hrs) 


3 

2 

4 

Sintering Temp. (*^0 


12^0 

1220 

1260 

Nb«0^ (atom fraction) 

2 5 

"3 

0.045 

0.01 

0.08 

MnO^ (atom fraction) 


0.025 

0.01 

0.04 

Al^O- (atom fraction) 


0.06 

0 

0.12 

Bi^Q^ ^3.tom fraction) 

^6 

0.01 

0 

0.02 


3* 3. 3 SELECTION OF EXPERIMENTAL DESIGN : 

Since there are six variables, a 2 fractional factorial 






design around the base levels was employed- The design matrix is 


shown in Table 4.1. The basis for choosing ^6 * 

as these combinations are expected to have minimum 
int er act ions - Experiments were carried out in a randomized 

sequence to avoid bias, on the part of experimen ter - The values 

of all the responses (dielectric constant, mechanical quality 
factor) for the first move of experiments (runs numbered 1 to 16) 
as planned in the design matrix are given in Table 4.1. 
Experiment No. 4 was repeated three times so as to obtain an 
estimate of the error. 

3-4 SAMPLE PREPARATION: 

PIT compositions with four dopants were prepB.ri3td by standard 
ceramic processing route. Dopants were Nb^O^, , Bi^O-^^ Al^O-^- 

The morphotropic phase boundary in Pb(Zr Ti )0^ lies B,t ^ 

0.535. Eight different compositions were prepared having this 
value of X but part of ZrO^, "^^^2 PbO were replaced by dopants 

so that composition can be represented by 

Pb, Bi {fzr- c:-r=Ti^ yi/cT A1 Mn Nb jO- 

1-m fn U 0*535 0. 465J l-y-z-6J y z ojJ 3 

Eight sets of samples of following composition were prepared. 


(1) 

[^^'^0.535 ^^0.465 

^0.79 

^^0.12 

""o.oi 

"^.OB ] °3 

(11) 

Pb ^^0.465 

^0.98 

’^‘^0.01 

"‘•o.oi 

]°3 

(111 ) 

[^^*^0.535 ^^0 

.465 

^0.85 ^^0.12 ‘^" 0.04 *^^0.01 ] °3 

( IV ) 

Pb 535 ‘'"^0.465^ 

0,88 

""0,04 

'^‘^0,08 

]°3 

(v) 

’^*^0.98 02 [^^’' 0.535 ^^0 

. 465*0. 

86 ^^0. 

12 ^'^O.Ol *^*^ 0 , 01 ] 



(vi) ®^0.02 


0.535 "^^0.465 ’o.91 ^^0.01 '^^O.OB ] °3 

(vii) ®^0.02 [^^'' 0.535 '*'’■0,465 ^0.95 *^'^0.04 *^*^0.01 ] °3 

<viii) PbQ^ 9 Q ®^o.02[**^''o. 535^^0. 465^0. 76 ^^0. 12*^"0.04*^*^0.08] °3 

Each of these above compositions were sintered at tv^o 
different temperatures ^i) 1220'^C and (ii's 1260'^C« The sintering 

time was 2 hrs or 4 hrs as per the design matrix (Table 4,1). A 
total of 16 distinct combination of all the factors plus three 
repeats of combination No, 4 were used giving a total of 19 
experiments as per Table 4.i« 5 wtVi excess PbO were added for 

those compositions which were sintered at high temperature 1260'^C« 
20 gms of batch materials for each of 16 compositions were 
prepared. Individual oxide component was weighed accurately. 
Then each batch were ball-milled for 10 hrs and dried in exactly 
same way as described in sections. 2. Dried powder was crushed in 
mortar pestle and then loosely poured in the recrystallised 
alumina crucible. Atmosphere powder was taken in a platinum 
crucible. It was placed over the loose powder. The system was 
covered with similar alumina crucible. The total assembly is 
shown in Fig. 3.1 (a). It was then put in the pit furnace. The 
calcination was done at 900'^C for 4 hours for all the batches. 

The calcined aggregate of each set was crushed into powder in 
agate mortar and pestle. The powder was then sieved through a 100 
mesh sieve and mixed with PVA ( Poly vinyl alcohol) solution (4*/# 
by wt) to yield 1 wt!4 PVA of the powders. Then water was added to 
make a slurry. After thorough mixing it was dried at 60^C ± in 

an oven. The dried powder formed a lump which was crushed again 



screen <ASTM)« 


in mortar pestle and passed through 100 mesh 
Pellets were JTiade fi?<ed amount of powder and well compacted 

using 5 tons pressure for 1 minute in a die of 12mm diameter* 


Pellets of 5 mm thickness were formed. 

The prepared pellets were sintered with atmosphere powder 
according to sintering schedule given in the Table 4.1. The 
heating rate was 3‘^C/min and the cooling rate was 4^C/min. The 
entire sintering assembly is shown in Fig.3.1<b)* Before 


sintering pellets were held at 600‘'^C for 4 hours to drive out 
adsorbed water and binder respectively. The weight of the pellets 
were taken after cooling down to room temperature. This is the 
weight of pellet before sintering. Finally the pellets were 
sintered and again weight of the pellets were taken after cooling 
down to room temperature- This is the weight of pellet after 
sintering. These two weights were required for finding out the 
amount of Pb - loss. These pellets were used for all further 
studies . 

3.5 XRD ANALYSIS; 

XRD analysis of sintered pellet for 16 different compositions 

were made by a REICH SEIFERT ISO DEBYEFLEY 2002 DIFRACTOHETER . 

CuK <X = 1.5405 A) radiation was used, 
a 

The X-ray diffraction plots of the sample were taken in the 
2.B range of 20"^ to 70*^ with the char’t speed 3 cm/min and scanning 
speed S^^'min and full scale intensity iOK counts/min. Also slow 
scanning were taken in the 2^ range of 30 —33 ^ 37 —39 , 43 —46 

With the following operating conditions. 



Atmospher powder 



Platinum crucible 

Recrystallised 
alumina crucible 


PIT powder 


Fig. 3. 1(a) Schematic diagrdm ot crucible arrange- 
ment during calcination. 



Alumina crucible 
AI2O3 powder packing 


Recrystallised alumina 

crucible 

Sample 

Alumina sheet 
Platinum disc 
Atmosphere powder 
Alumina powder 


Fig. 3. 1(b) Schematic diagram of crucible arrange- 
ment durina sinterina. 




Operating Conditions of Dif f r actojnaeter during X-ray Diffraction 


Current, Voltage 
Time constant 
Beam slit ui id th 
Detector slit width 
Scanning speed (2B) 
Chart speed 
Full Scale Intensity 

From the 
present were 
positions were 
3.6 


20mA, 30KV 
10 seconds 
1 «5 mm 
0.3 mm 

0.6 cm/min 
10 K coun ts/minute , 
20 K counts/minute, 
vs Intensity) 
peak positions. 
X-ray data files 


diffraction pattern <26? 
identified from the 
compared with standard 
MEASUREMENT OF SINTERED DENSTTY 


the phases 
The peak 


The sintered density of the prepared samples were measured by 
water displacement method. The initial dry weight of the samp 

was taken using a ACFOSET ER-120A ELECTRONIC BALANCE. The sample 
was then immersed in a water and kept under vacuum for 1 hour in 
order to allow water to enter into the pores of the sample. Then 
weight of sample suspended in water (w^> was taken. The sample w 
then taken out and dabbed with tissue paper and cotton. The 
weight of the sample, saturated with water was finally taken. The 
sintered density was calculated from the following relation. 

SINTERED DENSITY = -r, i-n- 

3.7 SCANNING ELECTRON MICROSCOPY s 

Fracture surface of the PZT samples were observed using JSM — 
840A Scanning Microscope, JEOL . 

The sample was broken before putting it in the microscope and 



was cleaned thoroughly using ultrasonic vibrator. Then it was 
mounted on an aluminium stub using a silver paint. For each sample 
exposed surface was coated with gold— pal ad lum in a dc sputtering 


^33 * 


unit (Anatech Ltd.) an operating voltage 5 K.V., plasma atomic 

discharge current 10 mA, and for 5-6 minutes. 

3.8 MEASUREMENT OF PIEZOELECTRIC PROPERTIES : 

3. 8.1 

Disc shaped samples of 1 mm thickness were cut using a 
diamond rimmed circular blade, polished with 3(0), 4/0 emery paper 
and cleaned with actone in an ultrasonic bath. The discs were 
then coated fay silver paint on both surfaces except the edges. 
Using a multimeter the continuity of each surface and the 
discontinuity between the two surfaces was checked. A strong 
electric field was then applied to the disc placed tightly between 
a brass cathode and an aluminium anode using a fixture 
schematically shown in Fig. 3. 2. 

To avoid dielectric breakdown, the whole fixture was put in a 
beaker containing silicone oil. The beaker was heated and 
temperature of the oil was maintained between 110 to 120*^0. The 
poling field was along the thickness. The fallowing conditions 
were used 

Applied field = 3.5 KV/mm 

Temperature = 110 - 120'^C 

Time = 30 mins. 

When the temperature reached about lOO'^C, the elect^'ic field was 
switched on and after holding at a temperature around llO^C to 
120*^0 for 30 minutes, the sample was cooled by switching off the 




Fig. 3.2 Schematic diagram of poling fixture. 



heater with the field on- When the temperature reached the roam 
value^ the field was switched off. Cooling to room temperature 
took about 7-B minutes. The discs were then taken out of the 
pol iHQ fixture and silicone oil was removed by Qently rubbinQ with 
cotton, care being taken not to remove the Ag coating. The poled 
samples were then kept for 24 hours after which measurements 

w w 

were made. A channel products, Belincourt d.^.^ meter was used for 

W 

d-- measurement, at 25 Hz. 

3. 8. S Measur ement of Resonance and Antiresonance Frequencies 
Mea s ur ement T ec hni ques 

Two types of techniques can be used: 

(i) Static method 
< 11 ) Dynamic method 


P lezoelectr 1 C 

constants 

may 

be 

measured 

by 

the 

techniques 

of 

static method with 

reasonab 1 e 

accuracy but 

the 

p rec ision 

IS 

inferior to that 

obtained 

with 

dynamic 

me thod . 

Hence 

the 


discussion is restricted to later method only. 

Elastic bodies show numerous resonances. But the most 
pronounced are those where the body can just accommodate one half 
wavelength of a standing wave. Piezoelectric effect is used to 
excite such elastic wave and to observe interaction of the 
mechanical ^resonance with electric field behaviour. 

Resonance Method 

This IS a simple dynamic method for evaluating piezoelectric 
constants by measuring the resonant (fr) and antiresonant <fa) 
frequencies corresponding to a longitudinal mode by means of a 
simple electric circuit as shown in the Fig. 3. 3. 




Fig. 3.3 Transmission net work for determination 
of fj- and fQ of piezoelectric ceramic 
vibrator. 



When the piezoelectric material is put in the circuit, the 
electric field sets up vibration in the material through direct 
piezoelectric effects The commonly encountered modes of motion of 
bars, plates and disks excited piezoelectrical ly by electric field 
are 

!• Length extensional mode 
11 * Overtone con tour-shear mode 
III, Contour extensional mode 
iv. Contour mode 

V, Contour extensional mode (disk) 

The first two modes of motion are for narrow bars, third and 
forth correspond to square plate and the last mode of motion is 
for disk. 

The periodic vibration of the crystal causes periodic 
piezoelect r 1 C changes on the electrode through converse effects. 
The properties of any mode of a lightly damped mechanical 
vibrating system (excited p i e zoe lec t r ical ly ) can be represented 
near resonance by an equivalent electric circuit as shown in 
Fig* 1.3, which consist of a capacitance , inductance , and 

resistance , in series, shunted fay a second capacitance in 

parallel. These are the fundamental parameters of the 
piezoelectric vibrator. The parameters are independent of 

frequency if the vibrator has no other mode of motion near 
resonance . 

The equivalent series resistance and reactance of 

R.L^C^C circuit is given by 
1 1 1 o 


the 





'? 'P '7 '7 

WC X + W'~C‘'R7 
o I Q 1 



r 2 2i 

X - WC 

Rf + X 

Q 

I i 

i - WC 

2 2 2 2 
X + W C^R 

0 

0 


where X = WL, - 


^ - WC~’ L^, 


branch in the 


series. The resonance frequencies are obtained by setting X = 0. 

At resonance R is nearly zero- 

Solving this series resonance, f s , is given by 
. ^ 1 


2?7(L^C^ ) 


1/2 


Similarly parallel resonance frequency can be obtained as 

1/2 

l t 

1 


C + C 


2n 


Q 1 j 

1-,C c, j 

^ 1 Q 1 J 


Variation of series reactance, X , with frequency is shown in Fig. 
1.3* Piezoelectric properties are measured using their resonant 
fr (minimum impedance) and an t i resonance fa (maximum impedance) 
frequencies. Complete scheme of measurements are given in 
Appendix 1. 

Here the measurements are restricted to disk shaped samples. 

Only fr and fa are measured by using 4t94A IMPEDANCE /GAIN PHASE 
ANALYSER (HEWLETT-PACKARD). 

Planar coupling coefficient kp , is obtained from the fr and 
fa values of the thin disk with faces perpend icular to z-axis as 


f ol lows ; 

. 2 


( i - ) J ^ ( 1 •+ -I- 


Af 


)3 - n , n + ~)3 

'1 fr o '1 fr 


2 

l-kp 


+ ^>3 


where J = Bessel function of fxrst kind and zero order 
o 

Jj = Bessel function of first kind and first order 


I.T.. 



'n 


= Lowest positive root of <l+<y^)J, = r;J (.n) 

^ l a 

For fy = 0.3 (poisson ratio), r/ = 2,05. 

The relationship is plotted in F 19 , 3.4 where kp is shown as 
a function of Af/fr C43. The value of kp is read directly frorn 
this plot, 
k. 
k: 


33 

and 

are calculated as 

follows : 

2 

n 

fr Cn AfV 

fl I 2 fTj 


33 

2 


2 

= V- 



31 

[ 

2 J p 


33 

was 

measured by using the 

CHANNEL PRODUCTS, 


BERLINCOURT PIEZO d^^ METER, 


IS calculated as follows; 


Si = *<31 


[^3 < 0 ^ 


/2 


coulombs /New ton 


where 


2 . 2 , 2 
rt d f r 




^11 '^<1 
d = diameter of disk 

p = density of disk 

= 2.05 

= 0.3 (usually between 0.29 to 0.33) 

T 

= Free dielectric constant, obtained as shown later 

Q--- IS given by 
33 . 

33 Volt meter 
9- 


'33 


T Newton 




(relative clamped dielectric constant, longitudinal) is 

33 

calculated as: 

S CxT 


K 


A X 8.854 X 10 


-12 


C = capacitances in Farad 
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Fig. 3.4 Planar coupling factor of thin disk. 



t = thickness in m 

'TV 

A = Area in 


9*854 X 10 Farad/meter 

s .^0 

T q 

(relative free dielectric constant) is calculated from by 

q 33 

T ^33 

^33 ** 2 

p 

Q(t\ can be determined by means of the piezoelectric resonance 


with a radial fundamental mode* Thus value of Qm is obtained by 
means of the relation given below, 

ffir = - (ti)''} 

where 


fr = resonance frequency 
fa = an 1 1 resonance frequency 
R - resistance at resonance 


C = Capacitance at IKHz 



CHAPTER 4 

RESULTS AND DISCUSSION 

4.1 DENSITY AND Pb LOSS : 

Theoretical density and */. Pb loss of the 16 sets of PZT 
compositions are shown in table 4,6. 

The results are plotted as a function of (a) different dopant 
concentration and <b) sintering temperature in Fig. 4.1. Here all 
the data points are used. This is a valid procedure in the 
analysis of data from the designed experiments. 

With the variation of doping concentration density varies. 
It IS observed that % theoretical density reaches higher values 
when doping with 0.08 atom fraction of Nb 0 in all but one 

jC, w 

experiment. Also with zero atom fraction of A1„0^, the density 
tends to reach maximum values. The effect of other dopants on 
densities is not discernible. Increasing the sintering 

temperature from 1220^C to 1260^^0 also does not result in any 
increase in density. 

The effect of dopant on the sintered density is through their 
effect on the rate of gram growth. If the grain growth is 

restricted then high density can be achieved because the diffusion 
distance for vacancies from a pore in the centre of a grain to the 
grain boundary is less. Thus the aluminium substitutes at B sites 
and creates an oxygen vacancy while Nb also substitutes at B site 
but creates Pb vacancies. Each is associated with its 
corresponding charge vacancies and seggregate to grain boundary 
making grain growth difficult Ct51. 

When both Nb and A1 are doped together, they associate with 



each other rather than with vacancies and so do not seQQreQate to 


grain boundary and do not affect gram growth* Individually Hb is 
found to be most effective and A1 less effective while when both 
are present no significant effect is observed* This is consistent 
with our observation that highest densities are obtained when only 
Nb IS present and the density is reduced when both are added 
together* 

The data for Pb loss has been similarly plotted in Fig*4*l* 

It is found that Pb loss increases with increase in 

sintering temperature from 1220'^C to 1260^C irrespective of the 
composition* Also variation of % Pb loss is significant with the 
variation of atom fraction of A1^0_ from O to 0.12* There is also 

jC o 

Pb gam for some composition* 

The effect of a trivalent dopant on Pb loss or gam has been 
discussed by Ulu et al* II343* The general equation for a trivalent 
dopant with some excess PbO (AO) is represented by the following 
equat ion , 

For trivalent doping: 


AB0« + XM-D- 4- 2pkA0 *4 ^ , 

3 23 < l+2pX~3qK- (3r /2) X 


2(q+r)x bJqx 2px *-^<r/2>x 3 2 


+(3q + XAO 


For the specific case of M on the A site with A-site 


vacancy we have 


ABO^ + '‘^(1-3X) ^2.x + 3 XaJ j loss 

For on A site with B site vacancy only 


For M on B site with 0 vacancy 







ABO + XM 0^ + 2XAQ (gain) -♦ CA 1Q_ O 

l+2x 2x 3+5x x/2 

Thus depending on the site to which the dopant goes and the 

site at which vacancy is created one can have either a loss or a 

gam in the weight- In general if the dopant goes to B site 

then an equal number of Pb atoms come to the newly create A site 

leading to an increasing weight- However even though it is known 
+3 

fchat A1 goes to the B site^ we find that there is a loss of 
weight instead of an increase in weight- Obviously the effect of 
other dopants dominate- 

4- a SCANNING ELECTRON MICROGRAPHS CSEM:> ANALYSIS : 

Hicrost ructure of six sets of samples are shown here (Fig- 
4-2-1, 4-2-2*, 4-2-3, 4-2*4)* From the microst ructures it is 
clear that samples (Ri, R2, R13, R14) which are poled of smaller 
in grain size and samples (R9, R12) that cannot be poled are of 
larger in grain size- 

The large gram size results in-excessive strain during 
poling causing the samples to break- Most of the samples (all 
except one) contain low amounts of Nb and also high amount of Mn . 
Apparently there is complex interaction between the various 
dopants which control the gram growth. The breakage of sample 
during poling appears to be solely because of large gram sizes- 
4-3 XRD ANALYSIS 

The obtained results (Table 4-6) indicate that the 
tetragonal ity increases on increasing sintering temperature for 
all compositions except one (R3, R5)* Some represent at i ve X— ray 
scans are shown in Fig- 4-3 and 4-4- Appreciable shifts m the 
peak positions are also observed for samples of different 
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composition <Fiq. 4.3). The other sample had 

naa Similar patterns and 

showed phases as listed in Table 4 . 6 . 

4.4 PIEZOELECTRIC PROPERTIES 

It has been observed that when the tetragonal ity decreases 
the average d 33 values (Table 4.7) get higher obviously because a 
greater amount of the sample mass has a chemical composition 
closer to that at the MPB. d 33 values are higher for the 

composition which does not contain AI 2 O 3 . Some samples broke 
during poling and could not be poled due to large strains 
generated during application of voltage. Dielectric constant 
values (Table 4.7) before poling are higher than after poling. 
Maximum dielectric constant obtained was 902 for R 6 . 

The variation of reactance with frequencies for the 16 sets 
of samples is obtained from impedance analyzer and some of them 
are shown in Fig. 4. 5, 4.6, 4.7 and 4.8. All the poled samples 
responded well to the impedance analyzer. The values of resonant, 
antiresonant frequencies and mechanical quality factor <Qm> are 


elastic compliance at constant electric 
electromechanical coupling coefficients (k^ 
1 ong 1 tud 1 fi a 1 voltage coefficient (g.^ 3 ), frequency 
are calculated and result are given in Table 4.7. 


range of 63 to 235. These values are low. Electromechanical 
coupling coefficient values obtained are in the range of 0.27 to 
0.45. 




constant 


field 


*^31 ’ 

**^33^ ’ 

constant 

(F.C. ) 

obtained 

are in 


In our experiment we have added four dopants and also their 



(3d) ^0 
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Variation of Reactances with frequencies 
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es with frequencies 


atom fractions are varied all at a timem There are 8 sets of 
different composition* There are no similarities between any two 
sets of composition* But amonQ the dopants some are donors and 
some are accepters* Thus Nb_0^ , Bi^O^ are donors and Al^O^ acts 
as an accepter* 

Donor doped material creates Pb vacancies* These vacancies 
facilitate domain wall motion with the application of very small 
electrical field or mechanical stresses* Thus it increases bath 

elastic and dielectric compliance and loss factor* Also vacancies 

+2 

in Pb position reduce the local stress in the domains* Thus 
domain can easily switch upon application or removal of the field* 
The cess of switching is related to relatively low coercive field 
and squareness of hysteresis loops- So mechanical quality factor 
Qm decreases with decrease of internal friction and hysteresis* 
Poling is easier for donor doped PZT material* As donor doped PZT 
materials show saturated normal hysteresis loops just after the 
application of field, they are ineffective in generating internal 
bias field and also ineffective in stabilising the pol ar i zat ion * 
In our experiment sets R2, R4, R8, R6, R12, RIO, R15, R13 are 
mainly donor doped PZT material* Thus we find that they have high 
dielectric constant ranging from 510 to 902 at IKHz* Values of Qm 
is also low except R2* This may be due to the presence of HnO^ 
which increases the Qm value somewhat by creating space charge 
which pins the domain wall motion and stabilizes the domains* Also 
Re, R12, RIO could not be poled* Poling is difficult due to the 
restricted domain wall motions* 

But in the case of accepter doped PZT materials vacancies in 



fche oxygen position will tend to shrink cell size. This shrinkage 
and distortion of cells is responsible for increasing mechanical 
quality factor Qm, increasing coercive force and lowering 
dielectric constant. Oxygen vacancies create large amount of 
space charges which slowly moves to the domain wails by diffusion, 
pin the walls and stabilize the domains. Thus domain wall motion 
IS suppressed and internal friction due to this motion decreases. 
Thus accepter doped PIT materials show unsaturated constricted 
hysteresis loops due to the existence of internal bias field. 
Poling IS carried on accompanied by the domain wall motion. So 
for accepter doped PZT internal bias field affects on the poling 
process . 

In the run numbers Ri, R7, R5 , R3, Rll, R9, R16, Ri4 both 
accepter and donors are present. They have low dielectric 
constant value. R5 , Rll, R9 could not be poled. But none of them 
has high mechanical quality factor (Qm) value. For filter 
ceramics mechanical losses should be as low as possible and 
consequently the Qm should be as high as possible. Mechanical and 
electrical losses in soft piezoceramics mainly arise due to 
vibration of domain walls and their interaction with lattice point 
defects. Mechanical losses in piezoelectr ic material is related 
to electrical conductivity ia) which is discussed in chapter 1-9. 

of Mn doped ceramics is far higher than that of ceramics doped 
with donors or complex oxides. Mn sometimes shows anomalously 
large dielectric losses due to inhomogeneous high conduct i v i t les 
at grain boundaries. So the value of Qm for Mn doped ceramics is 
low. Also Mn gives stabilising effect, i.e. it gives high 



stability with respect to frequency constant which makes them very 
useful in frequency filters. 

Electromechanical coupling coefficient values for all the 
sets are quite good. All the composition prepared are near the 
MPB. MPB IS of special interest because the permittivity and 
reach their peak values. For all PIT ceramics the dielectric, 
piezoelectric and elastic properties and their temperature 
coefficients are closely related to the various f erroelectric 
phase transitions. Spontaneous polarization in Rhombohedral phase 
is more favourable than in tetragonal phase. Also spontaneous 

strain is less in Rhombohedral than in tetragonal phase. 

4.5 OPTIMIZATION STUDIES 

MODEL FITTING 


The fallowing first order model has been fitted to the 
experimental data. 

Y = + b^x^ + ^ 3 X 3 + + b^x^ + b^x^ 

where b b, are the best fitted values of coefficients, x is 

a 6 

the ith variable in its coded form and V is the predicted value of 
the response. The coefficients can be estimated by the least 

square technique as 
T T 

b = (XX) X Y 



Table 4.1 


The 2 


6 -: 


First Order (First Movei 


Design Matrix 


and values of the Responses 


RUN 

NO. 

1 

y 

FACTORS 

IN 

CODED 

FORM 


VALUES OF 

RESPONSES 

1 1 

i 

1 

^2 

^3 

^4 

X 

<X 

X w 


DIELECTRIC 

CONSTANT 

(BEFORE 

POLING) 

MECHANICAL 

QUALITY 

FACTOR 

1 

i 

4 - 

+ 

4 

- 

4 

- 

623 i 

109.78 

2 

4 - 

4 - 

- 

- 

- 

- 

i 

522 

235.02 

3 

— 

- 


- 

4 

4 - 

579 

86 * 1 6 

4 

4 * 

- 

- 

- 

- 


510 

109. 13 

5 

- 

4. 


- 

4 

4 

405 

0 

6 

^ - 

- 

4 - 



4 - 

902 

63.40 

I 

7 


- 

4 

- 

4 - 

- 

598 

1 42.90 

j 

B 

- 

4- 

4 

- 

- 

4 

i 639 

t 

1 0 

1 

9 

- 

- 

- 

4 

+ 

- 

525 


10 

+ 

- 


4 

- 

4 

1 642 

1 

0 

11 

- 

4- 

- 

4 

4 

- 

451 

0 

12 

4. 

4 - 

- 

4 

“ 

4 

667 

0 

13 

- 

- 

4 

4 

- 

- 

540 

118,83 

14 

4 - 

- 

4 

4 

4 - 

4 

543 

117.02 

15 

i 

4. 

4 

4 

- 

- 

533 

81,59 

16 

! 

4- 

4 

+ 

4 

4 - 

Repeat 

4 - 

Trials 

502 

119.61 

17 

4 - 

- 

- 

- 


- 

486 

180,20 

IB 

4- 

- 

- 


- 

- 

520 

1 110.50 

1 

[ 

19 

4- 

- 

- 

- 

- 

1 

522 

109.60 

1 

1 


Aihere X is the design matrix. is the transpose of the design 

natriK. Y is the vector responses and b is the vector of 
roef f ) c 1 en ts . The details of the method are given by Draper and 
i^mith C361. This is then tested for model adequacy by a 
Lack— of— fit F test. With the experimental results as given in 
Table. 4-, 1, the fitted first order models for dielectric constant 
ind mechanical quality factor are given belotAi by equations (4.1) 
M^d (4.2>. 

^D.C. "" ~ “ 2.34x^ + 4.385X.J - 1 .58x - 3.79x^ + 


4.37x, 

o 

' = 7.92+3.e5x +0.57x.,+0.75x_-0.83x^-0.e3x^ + i .37x^ 

wm 1 2 3 4 5 6 


'5 

4.1 

4.2 


jhere represent the predicted values of dielectric 
•onstant, mechanical quality factor, respectively. The lack of 
it F tests for the above models are shown in Table 4.2 and 4.3 
“espec t 1 ve 1 y . A model calculation of the analysis of variance 
ANDVA) for conducting the Lack-of-fit F test has been done using 
LAG LIBRARY ROUTINE (as given in Appendix 2>. 

It IS observed from Table 4.2 and 4.3 that a simple 
irst-order model is not adequate for describing the dielectric 
onstant and mechanical quality factor data. Hence a higher order 
lodel incorporating cross-product terms have been proposed and 
ested for Lack-of-fit F-test. 

. 5 . 


This IS given in Table 4.4 and 



Table 4.2 


First order Response Surface Strategy (First Move) : 
Test for the Adequacy of the Model for Dielectric 


Constant 

A 

-D.C. "" - 0.56x^ - 2.34x2 - 1.58x^ - 3.79Xg + ‘^.37x^ 

ANQVA 


Source 

Degrees of freedom 

Sum Square 

Mean Square 1 

F- Value 

Model 

7 

24.5 

i 

3.51 I 

i 

4062 

Fit 

8 

0.409 

1 

0.051 1 

i 

59.1 

Error 

4 

0.0035 

1 

0.000864 

1 



Total 19 

] 

1 

1 

1 

1 



= 14.80 

Therefore, Model is not adequate. 

Table 4.3. : First order Response Surface Strategy (First Move) 


Test for the Adequacy of the Model for Qm 

Y„ = 7.92 + 3.85X, +0.57x„ +0.75><_ - 0.83x^- 0.83x-+ 1.37x^ 
Qm 1 2 3 4 5 6 


Source 

DOF 

Sum Square 1 

Mean Square 

F - Value 

Model 

7 

0.609 1 

0.087 

1000 

Fit 

8 

1 

0.126 

0.0158 

1 

182 

Error 

4 ! 

0.004 

0.000868 1 



Total 19 


Model IS not adequate. 

The cross product terms which have influence on dielectric 
constant value are given below with the values of the 
coefficients, 


b- = 572.5 b 

o 


-29.64 b.^ 


37.62 


- 22.01 


37.49 


35.12 






-44.14 


17.51 


-58.5 


Hence a better model can be proposed incorporating the above 
terms. The model then takes the following form, 

A 

^.C. ^ 572.5 - 29.64x^ + 37.62x., - 22.01 x. + 37.45 x, + 


Table 4.4 


+ 35,12 — 44.14XjX.^ + IT.Slx^x^ — SS.Sx^x^x^ 

First Order Response Surface Strategy (Second 
Move) : Test for the Adequacy of the Model for 


Dielectric Constant 



Therefore, model is adequate.. 

Similarly for Qm, proposed model even with cross-product 
terms does not satisfy the Lack-of-fit F-test. Anova Model with 
lowest F-value achieved is shown value, 

A 

Y_ = 78..55 + 8.11k^ - 7m65x. - 13.01 k, 30.27k.k-K- + 

Qm 3 4 6 1 4 5 

18 -28 k 21.19k^x,,x. 

4 6 2 4 6 

Table 4-5 : First Order Response Surface Strategy (Second Move) : 

Test for the Adequacy of the Model for Qm 


Source 

; OOF 1 

Sum Square 

1 Mean Square 

i F - Value 

Model 

7 1 

I 

69.1 

1 

9.367 

1140 

Fit ■ 

1 

8 ! 

6.48 

0.560 

64.5 

Error 

4 

- 0.0347 

0.00868 


Total 

19 





^ 0.01 


(4,8) 


14.80 








Therefore, model is not adequate. 

OPTIMIZATION? 

□f the above responses, dielectric constant of the sample was 
selected as the response to be optimized with respect to the 
independent variables. Thus data for the mechanical quality 
factor was also processed but no adequate model has been found 
using analysis of variance (ANOVA) . As a result it is not 
possible to propose any optimised set of independent variables, 
so that the response mechanical quality factor can be maKimised. 

Coded value of the independent variables of sample no. R6 has 
been fitted into the adequate model of dielectric constant. The 
value of response i.e. dielectric constant comes out to be 854.49. 
This value is very much close to the experimental value obtained 
(902> . 

To check the model further a sample was prepared taking all 

the independent variables as = 4 hrs, x^ = 12&0'^C, = O.Oi, 

X, =0.039, x_=0.12, X, = 0.02. The value of dielectric constant 
4 5 o 

obtained was 546. The result was also close to the value obtained 
from the model theoretically which was 559.94. 

Unconstrained maximization was not done for getting maximum 
value of dielectric constant. This needs to be studied later. 



Table 4.6 


% Theoretical Density, V. Pb Loss and Phases 


COMPOSITION SINT. SINT V. Pb V. PHASE 

DOPANTS CONCENTRATION TIME TEMP. (RUN LOSS THEORET. ANALYSIS 

(ATOM FRACTION) (hrs) NO) DENSITY 

Nb Mn A1 Bi (°C> 

0.08 0.01 0.12 0 4 1260 (Rl) 7. BO 96.83 R+T 

1220 (R7) 5.58 95.38 R 

0.01 0.01 004 1260 <R2) 4,36 91,86 T 

1220 (R4) -0.37 97.60 R 

0.01 0.01 0.12 0.02 2 1260 (RS) 5.52 88.95 R 

1220 (R3) 4.69 91.50 R+T 

0.08 0.01 0 0.02 2 1260 (R8) 2.71 98.11 T 

1220 (R6) -0.28 97.30 R 

0.01 0.04 0.12 O 2 1260 (Rll) 7.95 91.49 R 

1220 (R9) 2.46 94.10 R 

0.01 0,04 0 0.02 4 1260 (R12) 3.21 97.21 T 

1220 (RIO) -0,34 98.30 R 

0.08 0,04 002 1260 (R15) 3.03 96.63 R 

1220 (R13) -0.04 90.28 R 


0.08 0.04 0.12 0.02 4 1260 (R16) 7.82 94.97 R+T 

1220 (R14) 4.19 95.26 R 






ole 4.7 


PIEZOELECTRIC PROPERTIES 


SAMPLE 

DESCRIPTION 

AVG. 

d— — 

33 

(10“^^C/NJ 

DIELECTRIC CONST. 
<at 1 KHz) 

f r 

(Khz) 

fa 

(KHz ) 

Kp 

m 

BEFORE 

POLING 

AFTER 

POLING 

08011200/1260/4 

151 

623 

619 

242.6 

254.6 

0.35 

109.78 

08011200/1220/4 

138 

598 

436 

224,6 

235.4 

0.35 

42.90 

01010000/1260/4 

1 60 

522 

354 

197.0 

210.6 

0.40 

235.02 

01010000/1220/4 

179 

510 

279 

224.6 

244.0 

0.45 

109.13 

01011202/1260/2 

* 

405 

- 

- 

- 

- 

- 

01011202/1220/2 

123 

579 

464 

199.1 

204.4 

0.27 

86.16 

08010002/1260/2 

* 

639 

- 

1 - 

- 

• - 

- 

08010002/1220/2 

161 

902 

i 790 

1 

216.5 

228.5 


63.40 

01041200/1260/2 

1 

* i 

451 

- 


- 


- 

01041200/1220/2 

* 

525 i 


1 

... 1 

1 

” 1 


- 

01040002/1260/4 

* 

667 

i 

- 

- 


- 

01040002/1220/4 

* 

642 

- 

- 

_ 


- 

08040000/1260/2 

139 

533 

383 

242.0 

257.6 


81.59 

08040000/1220/2 

141 

540 

374 

248,0 

260.4 

0.35 

118.83 

08041202/1260/4 

142 

502 

373 

248,0 

260.3 

0.35 

119,61 

0804 1 202 / 1 220 / 4 

123 

543 

388 

250.4 

265.2 

0.36 

117.02 


Samples could not be poled. 
















Table 4.7 


Piezoelectric Properties (continued) 


RUN 

NO. 

33 


*<33 

*<31 

<10 C/N) 

<10“^^m^/N) 

F.C. 

<fr X r) 
KHz m 

R1 

705.41 

0.024 

0.34 

0.20 

53.15 

10.77 

1.191 

R7 

496,86 

0.031 

0,33 

0.20 

! 44.54 

10.74 

1.201 

R2 

421.43 

0.043 

1 0.38 

1 

0.23 

51.45 

12.90 

1.116 

R4 

349.84 

i 

0.057 

0.42 

1 0.26 

1 

45.20 

10.15 

1.222 

^3 1 

1 

1 

500.48! 

0.027 

! 

0.25 

1 

i 

0.16! 

1 

1 

37.87 ! 

12.95 

1.118 

16 

907.63 

0.020 1 

0.35 1 

0.21 ! 

63.81 

11.26 

1.161 

115 

440.03 

0.035 

0.37 

0.21 

42.74 

10.42 

1.212 

113 

426.21 

0.037 

0.33 

0.20 

41.40 

10.82 

1.226 

,16 

425.07 

0.037 

0.33 

0.20 

40.46 

10.36 

1.226 


14 


445.77 


0.031 


0.36 


0.21 


42.22 


10.04 


1.243 










CHAPTER 5 


CONCLUSION AND FURTHER SCOPE OF RESEARCH 

In the present experiments four dopants were selected for 
doping of PZT so as to study their effect with a view to arrive at 
compositions which may be suitable for application in filters. A 
design of experiments approach was used. The main parameter which 
was aimed to be optimized was the mechanical quality factor Qm. A 
set of 16 experiments were designed and conducted based on a 
fractional factorial design. 

It was found that out of 16 sets, m R5 , R8, R9, RIO, Rll, 
R12 sets, the grain size of the samples was very large so as to 
make it impossible to pole the samples. jpe Qm data for these 
sets therefore could not be obtained. For the purpose of 
analyzing the data using the response surface methodology (RSH) 
the value of Qm for the experiments was assumed to be zero. 
However, this did not help because the data could not be fitted to 
a model with reasonable degree of fit. Instead, the dielectric 
constant of unpoled samples was used as the variables to be 
analyzed using the RSM. This was successfully done. An arbitrary 
composition was then selected and a sample prepared and its 
dielectric constant measured. The measured value was found to be 
in reasonable agreement with the value predicted by the model. 
Thus it appears passible to design composition with desired 
dielectric constant using the model developed. 

This was an initial study to explore the possibility of using 
the design of experiments and RSH for PZT ceramics. Further work 
may be in the following directions! 



1 . 


The data should be used to carry 
maximization of the dielectric const 
value should be verified by experimen 
Another set of dopants should be sel 
one or two at a time and optimizat 
performed sequentially starting wi 
parameter followed by less important 


out and unconstrained 
ant and the predicted 
t - 

ected, preferably only 
ion studies should be 
th the most important 
parameters for filter- 
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APPENDIX 1 

Meftsureroeot .chetDe for pkroelectric 


ceramics 


Specimen and onenUtion 

iif ^ii» Bar plated on faces perpendicular to Z, 

Bar plated on face© perpendicular to Z. 

Bar plated on faces perpendicular to Z. 

Bar poled along Z, but plating on faces perpendicular to X. 


Heasurement scheme, or calculation method 


ft 


Sr plating on faces perpendicular to X. 

Bar poled along Z. but with plating on faces perpendicular to X. 

S^ Po|ed along Z. but with platmg c:. faces perpendicular to X. 
^ poled along Z, but with platmg on faces jieipendicular to X 
square plate with f^cea. perpendicular to Z. 


Square plate with feces perpendicular to Z. 

. d,,. Bar pla^ on faces p rpendicular to Z, and with length along Z. 


of the length 


Reeonaab© measurements 
cartonsional mode. 

From above, using «* (i _ ^ 

Measurement of low-frequency M^it^ce. 

Overtones of thickness aheax mode of the 

BMe bar as above, but with signal field 
along X. 

Measurement of low-frequency capacitance. 

Measurement of high-frequency (clamped) 
capacitance. 

= 1 — ei/ej. 

< = </(l - i„»). 

Fundamental resonance frequencies of the 
t^'o contour -extensional inodes of a 
square plate. 

*f. = -f, - iu**5. 

Beronance measurements of length exten- 
sional mode. 




Bar plated on faces perpendicular to Z, and with length along Z, 
Disk with faces perpendicular to Z, 

Disk with faces perpendicular to Z. 

Disk with faces perpendicular to Z. 

Disk with faces perpendicular to Z. 

Disk with faces perpendicular to Z. 

Calculation. 

Calculation. 


Resonance measurements of planar exten 
■ional mode. 

Besonanoo measurements of thickness ex- 
tensional mode. 

Measurement of response to hydrostatic 
pressure. 

== ^11* 

Overtones/^ of thickness extensional znode. 

Calculated from rf, and 


Calculation. 



1 



or 

Plate with faces perpendicular to Z. 


where 

^ 

Measurement of high-frequency (clam^^ed) 
capacitance. 


SloU: 8^^ #11 and d,i are generally negative, is taken to be positive when is negative. 
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APPENDIX - 2 

The program to analyse the data is given at the end of this 
Appendix. To run this program following is to be done; 

Imputing of’ data : 

The input data is to be entered in file sangital for 
dielectric constant and file sangita2 for Mechanical guality 
factor To input the data type the following command at the $ 

prompt . 


$ VI 


sangital i 

^ 4 — ! 


The screen will show the data presently stored in the file. The 
first few rows look like the fallowing: 


DATA FILE 

19 11 'U' 'M' 

+ 1 +1 +1 -1 +1 -1 +1 -1 623 

+ 1 +1 -1 -1 -1 -1 +1 +1 522 

-1 -1 -1 -1 +1 +1 +1 -i-i 579 
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Depending on whether the input file is sangital or sangita2 
the name has to be accordingly changed in the program. To access 
the program listing type the following: 

$ VI g02defe.f ^ 

Now move the cursor to the paragraph under "GET THE INPUT DATA" 
near the beginning of the program and change the file name as 
required. To come out of the listing, press 
ESc 

and then type, 

: wq J 

Mow to compile the program, type the following, 

t f77 +800 g02defe.f -Inag - 0 g02j 

'jow to run the program type the following, $ g02^ 

To see the output type the following, 
t Vi outl ^ (when input file is sangital) or 
& VI out2 ^ (when input file is sangita2) . 

The program does the first iteration using the constant SCO) 
ind the first coefficient B<2) provided in this example. In the 
lext iteration B(3) is also used and so on until the last 
iteration uses all the coefficients. The choice of coefficient 
zan then be changed if the F value obtained is not acceptable. 

4QTE : In the program the input value of dielectric constant 

and Q are divided by 50 for convenience, 
coefficient are therefore to be multiplied by 50. 


The 



vf *4 VI ^ ^ i ^ y * t’ ¥ h tr it ft i- f 

^ PROGRAM TO OBI AIN THE LINEAR LEAST SQUARE MODEL PITTING THE 
GIVEN EXPERIMENTAL DATA 
USES NAG SUBROUTINE GOaOEP 

PLEASE REPER TO DETAILED DOCUMENTATION IN NAG-14 REFERENCE MANUAL 
> 1 )>)>>> ^ >>>>>>>)>>>>>>>> X <<<<<<<<< <<<<<<<<<< <<<<<<<<<< < 

THE FOLLOWING IS A LIST OF THE DIFFERENT VARIABLES USED 


>>>>>;;> V >>> J >>>>>>>>>>)>><<<< < 

NMAX 

MMAX 

B(MMAX) 

COV 


X( NMAX I MMAX) 

1 i ^ *c * f 

implicit real ^8 
» .. Parameters 

INTEGER 
PARAMETER 
INTEGER 
PARAMETER 

♦ . . Local Scalars 
DOUBLE PRECISION 
INTEGER 
integer 
LOGICAL 
CHARACTER 
CHARACTER ^12 

♦ , Local Arrays . 
DOUBLE PRECISION 


<<<<<<<<<<<<<<< 

>> MAX NO. OF EXPTL. PTS. 

>> MAX NO. OF COEFFS. 

>> ARRAY OF COEFFICIENTS 
>> UPPER IRIANGLE OF THE COVARIANCE 
MATRIX 

>> ARRAY OF INDEPENDENT VARIABLES 

•» ^ I y 1 ^ f H 4 f ^ Ty ‘4 4 y y 4 H 4 444444 * hi 

( a-h , o-x ) 

MMAX, NMAX 
(MMAX=13,NMAX=E0) 

NIN, NOUT 
(NIN=25,N0UT=16) 


RSS, RSST, TOL 
I, IDF, IFAIL, 
IDX(MMAX-X ) 

SVD 

MEAN, WEIGHT 
INFIL 


INDX, IP, IRANK, J, M, N 


B(MMAX) , COV(MMAX^(MMAX+l )/£) , P ( MMAX-* ( MMAX + E ) ) 
+ Q(NMAX,MMAX+1 ) , SE(MMAX), UK ( MMAX^MMAX-^B^MMAX ) , 

UT(NMAX), XCNMAX, MMAX) , XX ( 3 ) 

^ . External Subroutines .. 

EXTERNAL F06FBF, GOEDDF, G0£DEF,FUN 

Executable Statements .. 

* GET THE INPUT DATA 


* WRITE (NOUT,^) 'MODEL FITTING RESULTS ' 

4 TYPE^' , ' QIME THE DATAFILE NAME' 

=» READ( N * ) INFIL 

^ URITEC=?^,*) 'THE INPUT FILE IS INFIL 

open<unit - nin, file = 'sangitaS') 
open(unit = nout , file ~ 'out£') 

4 4 ¥4 4 4 f 4 ¥4 y f ¥ ¥ ¥ f 4 ¥ 4 t 4 t 4 4 4 4 4 ¥ 4 4 4 4 4 4 4 4 4 4 ¥ ¥4444444 4 444^4*44444444^^^4* 


* Skip heading in data file 

44 4 ¥¥¥¥¥¥¥{ ¥ \ ¥ i 1 Kf«'» f 

READ ( NIN ; O 

READ (NIN,0 N, H, WEIGHT, MEAN 
=»- IF WEIGHT IS 'U', THEN UNWEIGHTED LEAST SQUARES DONE 


IF (N.LE NMAX .AND M LT.MMAX) THEN 

IF (WEIGHT. EQ. 'U' OR, WEIGHT. EQ 'w') THEN 
DO £0 I = 1 , N 

READ (NIN,^^) (X( I , J) , J=1 ,M) , Qd,!), WT(n 

4444,4 If ^ 4^ ¥ ¥ 4 4 4 4 ¥ * 4 ir ^ ^ ^ 4 * 4 4 ^ 4^ 4*4 *4 4* 4* 4 4 4 * 4 4 4^ 44 ¥ 4*** 4* ¥**** **** 

4 IN CASE OF MRR DATA, Q HAS TO BE DIVIDED BY £0 

44 ¥4 4 * ¥ * ¥ ¥ ¥ 4 ¥ ¥ ¥ ¥¥¥¥^'¥¥¥ ft" t 

* Q(l,1)=Q(I,n /£0 . OdO 

£0 CONTINUE 

ELSE 

DO 40 I 1 , N 

READ (NIN,’t) ( X ( I , J ) , J==1 , M) , 0(1,1) 

4 444 * r 4 4 ¥ 4 4 4 4 4 f 4' 4 4 4 54 4. 4 4 } ^ ¥ 4 4 4* 4 ‘4 4 4^ 4 ¥4 44 4 4 “» ¥ 44 44 4 ¥ 4 4 4 A- 4 4 4 4 4 4 4* 4 4 

4 IN CASE OF MRR DATA, 0 HAS TO BE DIVIDED BT £0 


QCI,I) = Q(I,1 )/50 OdO 

40 CONTINUE 

E^4D IF 

=4^ CALCULATION OF THE VARIOUS SUH SQUARES 
^ FIRST >> SUMSOUARES OF EXPTL DATA 

4C4..14 M-»^4^'4 4 4 ‘f4«4 r{444’44 4*4‘4‘4%4<^444.af:5Js4-4sHs|2Hii»c^+|*5f.=?*’4-^>»« + 4J«‘5|:%^%J^%s|!’!c^4jH.:«:j4.:4;a|6^ 

SSQH=0 . OdO 
do 1=1 , N 

SSQM-=SSOM + Q( I ; 1 ) ¥^2 
ENDDO 

^ UE OBTAIN THE MEAN FOR THE REPEATED RUNS IN THE PROCESS 

^ UE LOSE ONE DEGREE OF FREEDOM. 

4'l'44 4 44t>»-»4f41?-f4444>444444 44444H%4‘J»4-»|f44-54»f45444)K4'444'44.»4>»-4 4*%JH4. + 5|‘%>^34«»^%% 

SSQREP=0 . OdO 

3»^4«4. 4'4‘4'ftt444*4'f5r5 4444 J4444414v^4544f4 4’4t4.4%4‘4.4!4H ».4c3f4>K4:4+%3f.sj:3^5i!)4:%4e4cs|e3t:i|e34!»4!^ 

^ IBEG IS THE STARTING RUN NO FOR THE REPEATED RUNS 

34*4 4 4' 4 4 4*f4'4 4 4 ) ^ -> 4 4 4 *44 ^ ♦ 4*4 4'4'4-4‘t4-i424-4-4 %4«‘4’5|'4*4i^4^44 4+4-4 4+^4: + + + i4-3|;+*+- + + a|ej>;s»BJ|:4t 

IBEG •= 15 

I END = N 

do L = IBEG>IE^4D 

SSQREP = SSQREP + Q( I , 1 ) 

ENDDO 

jf4:4t4*4'+ + 4'4 4 4 4 4 4 fr 4 ? 4 + 4 4 4 4 4+4 + 4* 444 +4* + + 4 + + + 44 + 4 4 + + ++ + 4* + + + + + + + » + + + + + + + ♦♦ 

» THIS IS THE MEAN 

+ + + 44 + 4 + 4 + 4 + t + ’»< + 4 + 4 4 4 ++4 + + + 4 + ++ + + + + + + + + + 4 + + + + + ++ + + 4 + + + + + + + + + + + + + + + + 

SMEAN = SSQREP/f loat C IEND-IBEG+1 ) 

4^ 44:|,»^P444.4fi»^*f44+44+ + + + + + + + 4'4+4*++t + + + 4’ + + 4+ + + + ++ + + + + + + 

+ NEXT UE EVALUATE THE SUMSQ . OF REPEATED RUNS 

+ + + 4 ^44^4+1+-5 + ! + ^4 4 W+ + 4+ '>+44+4‘ + + + + 4‘ + + + 4 +444+ + 4++4 14- + 4+ + + + + + + + + + + + + 

SSQREP=0 . OdO 
do i=IBEG,IEND 

SSQREP = SSQREP + < Q ( I , 1 ) -SMEAN ) ♦♦E 
E^4DDO 

+ + t4 4 4:4^f4+-3 4 4' + ++ 4++!f f+ + + 4 M J.4+ + + + + + + + + + + + + 4 + + + 4 + + + + + + + + + + + + + + 4+ + + + + 

+ THIS IS THE EXPERIMENTAL ERROR DEGREE OF FREEDOM 

+ +4-4++444144'f4« + J44 4>l*4444+4’444+ + 4+ + 4+ + + + + + 4 +4++ + + 4+ + 4>+ + + + 4:+ + + + + + + 

DOFERR = f loat ( I END- 1 BEG ) 

!>::*: + !?+ + 4* + + '} 44 44 i+j + 44 + + 4+ + + + +^444+4+ + + + ++ + + 44+4* + + + 44+ + + + + + + + + + + + + + + + 

♦ THIS IS THE TOTAL DEGREE OF FREEDOM 

+ 4+ 4»+4|4 44 I * * 'i 4 *44 ^44 + 4+ 4 + 44 ^^^5j:t^4 4 4+ + + 4+ + +f + + + + + + +♦ + + +♦ 

DOFTOT=FLOATCN) 

4 + f‘4 4 4 4 f 4 4 4 ^4-4.4;4M 4 + 1 4+’4i>: + j^ + +3^+3r + + + + + l54+ + 4* + 4*4+ + ++'+4 4+ + + + + + + + + + + + + + + + 

+ UE CALCULATE THE MEAN SQUARES FOR THE TOTAL AND THE EXP. ERR. 

+ + 44t+4t 1441 1 + 7 I ♦»4’i444 4’+444*f*444+++444+ + ’»H’«- + 4+ 4+ 4 4-’4jf;»f:%4i*:4:^t)*,4.>fc + + >tcHji|e3H 

AMSQr=SSQTOT/DOFTOT 

AMSQE=SSQREP/DOFERR 

j|:3r4 4 + 4 4 + f 4 4 t 4 4 > 4 4 4 + 4 4 4 '4 4 +4'4+ + +44+3|4 +4H 4 4 1'4 +4 4 + + ++ %4 +4C + + + + + + + + + + + + + + 

+ S o t t o 1 1 ? r a in c e 

3|t4. 4-4t + 44 + 1+ 44 44 ^ 4 ‘4 4 4 4 4444+44+4+ + + + + + + +4 + + 4+ + + + + + + + +*+ + + + + + + 

TOL = 0 000001D0 
IP = 0 

IF aiEAN EQ 'M' .OR, MEAN EQ.'mM THEN 

+ 

CALL FOSFBFiN, 1 . 0D0,X( 1 ,MMAX> , 1 ) 

+ 

IFAIL = 0 

+ 

CALL GO£DEF(UEIGHT,N, I P , Q , NMAX , P , UT , X C 1 ,MMAX) ,RSS,TOL, IFAIL) 

* 

IP = 1 
END IF 
IDXC n = 0 
II = S 


=♦«+>> ^ ^ ^ t ? 4 4 if 5 4 5 -f f It \ ♦ j fif l^’f43f44H5»«4-4^‘4 4 4-f444’^‘'^44f + 4»4i>>»:344.a4if:j|;i>:t!He4e 

* UE READ THE COEFF TO BE ADDED 

60 READ ( ^4IN, • ) IIMOX 

IF (INDX G1 0) THEN 
IFAIL = -1 
ID.K(rn -r lUDX 

ic 

CALL G0£DEF<UE1GHT;N, I P , Q , NNAX , P , WT , X ( 1 . INDX) ,RSS,TOL, IFAIL) 

lie 

IF ( IFAIL EQ 0 ) THEN 

UPITE<nout,99993) 'CASE ' , IP 
IP = IP + 1 

WRI IE (NOUT , 

WRITE (NOUT, 99999) 'Variable', INDX, ' added' 

RSST = 0 ODD 
IFAIL = 0 

* 

CALL G0£DDF(N, I P , Q , NMAX , RSST , I DF , B , SE , CO V , SVD , IRANK,P, 

+ TOL,WK, IFAIL) 

« 

IF <SVD) THEN 

WRITE (NOUT,*) 'Model not of full rank' 

WRITE (NOUT,^) 

END IF 

f ¥ Hf f ¥ Jt4tt41'-44)44^-|-4't1et5|c3|e3|e‘1esfc54j*tf4 t t 4 '?«+:i‘5«£=|e4c|eHc'4:|<:>|ti« 

* UE OBTAIN THE RESIDUAL SUM OF SQUARES 

3|t:^{:4 ¥ ^ 4^ ¥ 4f ¥ 4t ¥ ir ¥ ► f f t ¥ ¥ ¥ f ¥ ¥ 4 4r f f ff1rici{5fc^J4*-4441f le-t 4 ¥'¥'¥* *’ |t?|r'|tie1e5r"f‘-4s|c*f:ie*i* 

WRITE (NOUT, 99998) 'Residual sum of squares = RSST 
WRITE (NOUT, 99999 ) 'Degrees of freedom idf 

WRITE (NOUT,’«‘) 

WRITE CNOUT,’^^ 

+ 'Parameter estimate Std, error t-value' 

* THE DEGREES OF FREEDOM FOR THE MODEL = r40 . OF COEFFICIENTS 

4^ 4t 4 ¥44 ¥4 i f t-l* ¥ ¥ ¥4 fi t- f ¥ f- ¥ ¥4 ¥444 ¥4 ¥4444 |c-*4< ¥ ¥ ¥ ¥ f. -ft 

D0FM0D=FL0AT< IP) 

WRITE (NOUT,*) 

DO 80 J “ t , IP 

•is** * * ¥444 ¥ A ^ 4 ¥'i 4 4 ¥ f 4 ¥44 ¥4444444 4 44444444*4 

4 HERE WE OBTAIN THE T- VALUES 

444 4 ¥ ¥ * ¥ |tf ¥ ¥ ¥ ¥ ¥ ¥ ¥ 4 ¥ ¥4 ¥44 ¥¥4444¥44 4 4 ¥4¥¥4 4 4 444¥ 4 44 44¥44 444444 

TTEST = B(J)/SE(J) 

WRITE (NOUT, 99997) IDX(J>, B(J), SE(J), TTEST 
80 CONTINUE 

4444¥444¥A¥^¥ » « 4 4 44 1ft4-5r’rH4f!444«4!4«4.>|^t.4t3lt4its!4:H«4c54'4:’|.i|!*li 

* WE CALCULATE THE REMAINING SUM SQUARES 

* 1 e, SUMSQ OF MODEL 

nt^ ^ 4 4 4 4 4 4 4 4 4f444iiii4:¥‘:¥44 4 44444^ 4 4 4' ^4 4 T4*4 4 ^44444*4 4 4444444 

SSQMOD«»SSQM-“RSSr 

4t45'V4444tf»4444 M4'4it1t4M46*44^-l*»4'4)4>»‘44J4t4«»l<4v4 I 4 4* 4 444 4' 41 >-4 4jf4«4454f454<4.»4J4i’l4>l«i|i 

* ^ SUMSO. OF FIT (Br SUBTRACTION) 

4^ 4 ¥ 4^ ^ 4c4 tl-4 4 4^ 4* 4 f 4c4 4 4-4 -»-r44 ¥4 44 444^ 4 ¥ 44 444 4H 4-4 if 444 4 4 4^4 4 44 4 4444 4*4444 

SSQLFT=RSST-SSQREP 

DOFRES‘=DOFTOT"DOFMOD 

DOFLFT«DOFRES-“DOFERR 

si£4cf4 4 ’4 4 ^ 44,4 4 ¥ V4- \, ¥ A ^ *1. 4 * *54 1-*4i**i*^4 4* 4-** 4^ *%******** 

NOU UE TACKLE THE MEAN SQUARES. 

i(.4!j44.4.4'}-4.4:J4444 444:4 4'4. 

AHSQH^SSQHOD/DOFriOD 
AMSQR^RSST/DOFRES 
AMSQL-SSQLFT/DOFLFT 
AMSQT^SSQM/DOFTOT 
FMOO = AMSQM/AMSQE 

3|*f4 4’444 \ i4444t‘i4 44 M4 i •44444*4 4’4^4c + 5|4lt4*4 444f44443*4%44*4:44.4 4'*4.***** 

* UE OBTAIN THE F- VALUE 


4 + 4 ■f‘»4 ^-4 *1 }<444444444.44%444 ■4'f’*4'»'l44'f3|1i4 l<4‘44.4‘4-4344‘4i’*-4-4‘ 

FVAL= AHSQL/AMSQE 

4 44 5V4 4 « 4^4*4 44 4 4 I 4 ^4 4 444 4 4 4444+44-4f44 4.*4'4 4 4 4 44c4^4t44:44444.»):4<4«44i*4<5»"J|c4!^* 

>!• AND COfJSTRUCT THE ANOVA TABLE 

!(|«4 4‘4*4tA4‘4 f444 4M‘4 4>4 4’4 4'4'4l4«*4'4444 4s4c453j.4s4f:j.:fa|}>-^4^44.,j.^5(i|,if^rf'^^^ 4.4c44:4t4:4:4e4<1c 

URI recnout .O' 

Write(nOUt |99994)' Ie4ef4fl'4'|r4t-ie4|cvf444f4f4t4t4«4:45+4^it4<4«*' 

URI TE(nout . 99994) ' ANOVA FOR THE PRESENT MODEL' 

Write(nOUt. 99994)' ‘4t4i4^4«r:{t|i3f.3f:4t‘f4l't4c4«4i5fts|c4t.4e*4t*4«*HI'ie' 

writeCnout. 99994) 

4 c 'Source DOF Sum Square Mean Square F-Value' 

wr ite(nout ; 99995) 'MODEL'. DOFMOD. SSQMOD. AMSQM , FMOD 

wrxteCnout . 99995) 'FIT DOFLFT. SSQLFT. AMSQL . FVAL 

wriTe(nout , 99996) 'ERROR'. DOFERR. SSQREP. AMSQE 

write(nout, 99994) ' 

•I' ' 

wnteCnout . 99996 ) 'TOTAL', DOFTOT, SSQM, AHSQT 

write(nout , 9 9 994 ) ' 

4 ' 

URITE( nout , 99994) ' 

ELSE IF (IFAIL.E0.3) THEN 

WRITE (NOUTO) ' ♦ variable not added *' 

ELSE 

GO TO 100 
END IF 
II « II + 1 

>|:4:^:3Jc 4::|;4:4:4:4!4!4;4’4:4 1:4. 4*4 t 4 ^'4. 4:4-4;4 H: ^ 4. ^ 

* PRESENT CASE OVER; NOU ADD THE NEXT COEFFICIENT AND START AFRESH 

GO TO 60 
Er4D IF 
END IF 

100 CONTINUE 

ft Vf ^ t f'4f4t4« + + :<v*iet:«c4‘'<'faf-4t3«cst:’f4!ie45^ic:{<ic 1c3f-4^4c4c^(4!4c:(c:|c3((9|; 

* IF YOU WANT THE PREDICTED VALUES OF THE DEPENDENT VARIABLE 
4C THEN THIS IS NEEDED 

:4c 54c + t4:BE C AR EF U L 4^ ic 4e sf: 

ic MODIFY THE SUBROUTINE FUN ACCORDINGLY BEFORE PROCEEDING WITH 
4C THE NEXT PART OF THE PROGRAM 

NUMPT ~ 30 

typefr.'DO YOU WANT TO OBTAIN THE FUNCTION VALUES?' 

Type^, ' 1 >> YES; B>> NO' 

r ead i ^ I RESP 

IF( IRESP EQ . 1 )THen 

^4't'f4'*4':44t 44 444 4 44 444>l44 4*4 4:'4i543f'4«4 3|'4'44e'4-4“4^t444c4-4 4'4'4-5*.^‘5?4?4c:*:*:4'‘4e4c4:4:5t:4c^»* 

THE RESULTS ARE STORED IN A FILE CALLED varl 
open(unit=£6,fi le='var1 ' ) 

101 type=r, 'UHICH VARIABLE YOU WANT TO VARY '7 1, £, 3, 4, 5, 6?' 

read ( ♦ , * ) I VAR 
if(3var eq 1) THEN 

Type-K . 'WHAT ARE THE VALUES AT WHICH YOU FIX X(£), X(3>, X(4), 

1- y(E) , X(6) ' 

RSAD(<','>) XX(£), XX( 3 ) , XX ( A ) , XX ( S » , XX C 6 ) 

t I S5X(£) 

t S =• X X ( 3 ) 
t3 = xx(4J 
t 4 = X X ( 5 ) 

T 5 = X X ( 6 ) 

ELSE IP (IVAR .EQ £) THEN 

typet.'WHAT ARE THE VALUES AT WHICH YOU FIX X(1), X{3), X(4), 
+ X ( 5 i / X ( 6 ) ' 

PEADCt.x) XX(1), XX( 3) ,XXC4 ) , XX(5 ) ,XX(^ 
t 1 = X < ( 1 ) 

t £ =■ X X ( 3 ) 


t 3 “ y X ( 4 ) 
t4 XX (51 
tS xx(6) 

ELSE IF (IVAR EQ 3) THEN 

type^.'UHAT ARE THE VALUES AT WHICH YOU FIX X(n, X(£), X(4), 
+ X(5) , X(6) ^ 

READ ( 41, IK ) XX( 1 ) , XX(2:) ,XX(4) , XX(S) , XXC 6 ) 
t 1 = XX ( 1) 

t£ = XX ( S ) 

t 3 == XX ( 4 ) 
t4 = xx(5) 
t 5 “ X X ( 6 ) 

ELSE IF (IVAR . EQ . 4) THEN 

tYpe4c,'UHAT ARE THE VALUES AT WHICH YOU FIX XCI), X(£), X(3), 
+ X( 5) , X( 6 ) ' 

PEAD(^,^) XX(1), XX(2:).XX(3),XX(5) AX(6) 

t t ^ XX ( I ) 
t £ •= X < ( £ ) 

t3 = xx(3) 
t4 = XX ( 5 ) 
t 5 = X X ( 6 ) 

ELSE IF (IVAR . EQ . 5) THEN 

typet^'WHAT ARE THE VALUES AT WHICH YOU FIX X(T>, X(2), X(3), 

+ X ( 4 ) . X ( 6 ) ^ 

READ( 4,0 XX( 1 ) , XX(£) ,XX(3) ,XX(4),XX(S) 

t t = X X ( 1 ) 
t £ = X X ( 2 ) 
t 3 = X X ( 3 ) 

t4 - XX ( 4 ) 
tS ==• xx(6) 

ELSE IF (IVAR EQ G) THEN 

type'40WHAT ARE THE VALUES AT WHICH YOU FIX XTI), X(2), X(3)| 
+ X ( 4 ) , X ( 5 ) ' 

READC^o) XX(1), XX(2).XX(3) ,XX(4) ,XX(5) 

t 1 = X X ( 1 ) 

t£ “ xx(£) ^ 

t3 == xx(3) 

1 4 - X X ( 4 ) 
t S » X X ( 5 ) 

END IF 

5ft 

t , ^VARYING xOIVAR,' FROM E.OdO To 4 OdO' 

v;r i te( 26 ; 'fr iNUMPT+1 
write(26, 99991 )t OtS 
do k k 1 = O NUMPT+1 

xX(ivar) = £ OdO-fS ODO=f DBLE(KK1 -1 ) /DBLE ( NUMPT ) 
c 1 1 f un( y , Xx , MMAX , B , 6 ) 
wr j te(2G.9999£)Xx< ivar),y 
ENDDO 

typefOWANT TO REPEAT THE FUNCTION CALCULATION?' 
t vpeH O 1 > > YES. £> > NO' 
r e aD ( 4 ,4 ) I RESP 
IP ( IRE5P EQ 1 ) THEN 
GO TO 1 0 I 
ENDIF 
END IF 

STOP 

99999 FORMAT (1X,A04.A) 

99998 FORMAT MX,A,D13 4) 

99997 FORMAT ( 1 X . ' B ( ' , I 3 ; M O 3D 1 5 4) 

99996 FORMAT (1X.A.4X.F3 0.£C£X,D1£ 3)> 

99995 FORMAT (}X,A.4X.F3 0 , 1 ( £X , D 1 2 . 3 ) ) 

99994 FORMAT M X , A ) 

99993 FORMAT ( 1 X . A, 121 

99992 FORMAT (1 x , F6 3. 5x , DIE 6) 



99991 FORMAT ( 1 x , f 4 1,1X,f4 1) 

99990 FORMAT (1X,A,D15.4) 

END 

♦ f f i ), t «■ f f )r f.f (tt* f-l< + Jet f+t'-ft-Jt Jtii t Jc-fit(c:»-t-'t--J= 

subroutine FUN ( y f un , x , MMAX , B£ / NVAR > 
implicit real *8 Ca-h,©-!) 

DIMENSION Ba(MMAX) ,x(NVAR) 

CSS + 

rFUN = BS(1) 1- BSCSl-rXil) + B£(3)>i'x(4) 

+ B£(4)4-x(5) + BS(5)»X{6) 

+ + B£( 6 ) ♦•X( 1 )>hX{S ) + BS(7)>»X( 1 )*X(4)H-X(5) 

“*= + B£(8) ♦'X(4)+X(6) 

* YFUN = esd) + B£C£)>kX( 1) + B£C3)<-x(3) 

* ■■¥ + BS(4) + x( 1 )*«•£ + B£( 5 I’^XI £)>!•*£ 

» % + 82(6 I’J'Xt 1 )**3 + B2(7)>f-X(3)»=*3 + B£ ( 8 ) + x d ) +*4 

RETURN 
END 

t. 4 + 4- 4' + 4 4 J J-**^^**** 4 4-4-4!4*‘4;4.4:!t:!j!4 4c4si;'t- + =“44 4 4:44»***4*** 



CASE I 


Variable 5: added 
Residual sum of squares = 
Degrees of freedom ~ 17 


Paramet er 

B( 0) 

B( 2) 


est i mat e 
1 1£££“i0£ 

- 3607E+00 


781 0Et0£ 

St d error 

. 4980E+00 
. 4980E + 00 


t -va 1 ue 

.a2B2E+02 
- . 7£43E-i"00 


ANOVA FOR THE PRESENT HODEL 




Source 

MODEL 

FIT 

ERROR 

DOF 

£ 

13 

4 . 

Sum Square 
. 24£Et04 
. 778Et02 
346E+00 

Mean Square 
. 121E+04 
.598E+01 
.864E-01 

F-Va 1 ue 
. 140E+05 
.692E+02 

TOTAL 

1 9 

,249E+04 

. 1 31 E+03 


CASE £ 

Var X ab I e 
Res i dua 1 

3 

sum 

added 

of squares = 

6446E-^02 



Degrees of freedom = 16 


Par amet e r 

est 1 mat e 

Std error 

t - val ue 

B( 

0 ) 

1 133E4-02 

4707E+00 

2408Et0£ 

B( 

H) 

~ . 47S8E+00 

4707Et00 

- . 1 01 7E+01 

BC 

3) 

.S662E+00 

.4707E+00 

. 1840E+01 


if ^ it: if fit 

ANOVA FOR THE PRESENT MODEL 
•»!* + ■(!.» ♦*********♦*!)! + 


Source 

DOF 

Sum Square 

Mean Square 

F*-Value 

MODEL 

3 

,243E+04 

, 81 OE+03 

. 938Et04 

FIT 

12 

. 641E + 02 

. 534E+01 

. 618E + 02 

ERROR 

4 

. 346E+00 

. 864E-01 


TOTAL 


249E+04 

. 131E-t03 



CASE 3 





Var iab le 

4 

added 



Residual 

sum 

of squares = 

. 6294E+02 


Degree s 

of freedom = 15 



Parameter 

est i mat e 

Std error 

t - V a 1 u e 

e< 0 > 


1130E^02 

4839E-+-00 

.2335E+02 

Bf 2) 

- 

.4438E+00 

4839E4-00 

.9172E+00 

B( 3) 


.901EE+00 

, 4839E+00 

. 1 862E+0t 

Bt 4) 


.£913Et00 

,4839E-J“00 

. 6021E+00 

¥ f 

f A ii * 



Ar40VA 

FOR THE PRESENT 

MODEL 


1 1 f -f 

i t t f t '-f if -ft ’ft #: f It tft t ff sf 


Source 

DOF 

-Sum Square 

Mean Square 

F-Value 

MODEL 

4 

, 243E +*04 

, 60SEt03 

.704E+04 

FIT 

1 1 

. 6e6E+oe 

.S69E+01 

, 6592+02 

ERROR 

4 . 

346EtO0 

.864E-01 




TOTAL 


, £49E+04 


. 1 31 £■^03 


1 9 


CA5<E 4 


V Hi r i a b 1 e 

c 

n <” j y •< .r* ri 



Re 

siduai s 

urn 

of squares « 

.S1£7E+0£ 


Degrees of 

freedom “ 14 



Paramet er 


€ St 1 mat e 

Std. error 

t -value 

B( 

0 ) 


. 1 1 39E’f 0£ 

. 4547E + 00 

.£504E+0£ 

ec 

£) 


5308E+00 

4547E+00 

-.11 67E+01 

BC 

3) 


, 81 4£E4-00 

4547E+00 

. 1 791E+01 

BC 

4) 

- 

3783E+00 

4547E+00 

-.83£0E-»*00 

er 

6 ) 


81 1 7E+00 

4547EH-00 

1 785E+01 




ANOVA FOR THE PRESENT MODEL 


Source 

DOF 

Sum Square 

Mean Square 

F-Value 

MODEL 

5 

£44£+04 

489E+03 

5E6E-^04 

FIT 

1 0 

. 509E+02 

509E+01 

589E+0£ 

ERROR 

4 

346E+00 

.8S4E-01 


TOTAL 

~Ti' 

£49E+04 

1 31 E-f03 



CASE 5 

Variable 7 added 

Residual sum of squares = 4275E-^-0£ 

Degrees of freedom = 13 

Parameter estimate Std. error t-value 


B( 

0) 


. 1 14EE+0£ 

.43£9E4*00 

.£846E+0£ 

Bi 

£) 

- 

S98£E+00 

.4389E+00 

. 13e£E+0t 

B{ 

3) 


.7468E*f00 

. 43£9E+00 

. 172SE+01 


4) 

- 

.44E7£'faO 

.43£9E+00 

1 030E+0t 

BC 

6) 


.7443E'^00 

.4329E+00 

. 1719E+01 

BC 

7) 


6968E+00 

.43£9E+00 

. 161 OE+01 







ANOVA 

FOR THE PRESENT 

MODEL 




It t -(c (c •* t 



Sour ce 

OOF 

Sum Square 

Mean Square 

F-Value 

MODEL 

6 

245E+04 

409E+03 

. 473E-I-04 

FIT 


9 

484E+08 

, 471 E + Ot 

. 545E+02 

ERROR 

4 

346E*«-00 

.864E-01 


TOTAL 

iT"* 

£49E*^04 

. 1 31 E+03 



CASE 6 


Variable 8 added 

Residual sum of squares = .3115E+02 

Degrees of freedom = 1£ 


Paramet er 

estimate 

Std error 

t -value 

BC 

0 ) 

t 1 38E+0E 

.3e61E+00 

.£948E'^0£ 

BC 

£) 

- S£6£E400 

.3861E-^00 

- . 1 363E+01 

B( 

3) 

. 81 8SE+00 

.3S61£t00 

.£1£1£+01 

BC 

4 ) 

~ 3T37E+00 

. 3861 E+00 

- . 9679E^‘00 



B( 

6) 

8163E+00 * 

.3e61E+00 

.£1 1 4E+01 

B( 

7) 

7688E+00 

.3861E+00 

. 1 991E-4-0t 

B( 

8) 

- 8168E+00 

. 3861 E+00 

-.£1 14E+Ot 


f ii f is if 

ANOVA FOR THE PRESENT MODEL 


Source 

DOF 

Sum Square 

Mean Square 

F-Value 

MODEL 

7 

E46E+04 

. 35£E-f-03 

.407E+04 

FIT 

8 . 

. 30SE+02 

. 385E + 0t 

,446E+0£ 

ERROR 

4 , 

346E+00 

,864E-01 


TOTAL 

Ts 

249E+04 

. 1 31 E+03 



CASE 7 


Variable 

1 0 

added 



Res 1 dua 1 

sum 

of squares = 

899SE+02 


Degrees 

of freedom = 1 1 



Faramet er 

est 1 mate 

Std. error 

t -va 1 ue 

B( 

0 ) 


1136E-^0£ 

. 3968E+00 

.£883E-f0a 

BC 


- 

5049E+00 

3968E+00 

. 1£73E+01 

B( 

3) 


.a401E+00 

. 3968E+00 

.£1 17E+01 

B( 

4 ) 

- 

.35£4EH'00 

.3968E+00 

. SBSeEi^OO 

B( 

& ) 


8376E4.00 

3968E+00 

.£111 E-f-OI 

B( 

7) 


r901E+00 

. 3968E + 00 

t991E+01 

B( 

8) 

- 

,7949£+00 

.39G8E+00 

.S003E+01 

B( 

1 0 ) 


£6£4E+00 

.3988E+00 

. 681 4E-f00 





ANOVA 

FOR THE PRESENT 

MODEL 



+ f- t' f- f ie t iff t if aft H t 4^ * 


Source 

DOF 

Sum Square 

Mean Square 

F-Value 

MODEL 

8. 

.£46E^04 

308E+03 

.357E+04 

FIT 


7 

. £96E+02 

4E3E+01 

.490E+02 

ERROR 

4 . 

346E+00 

.864E-0t 


TOTAL 

1 9 

£49Et-04 

. 131E+03 



CASE 8 





Var 1 ab le 

1 1 

added 



Residual 

sum 

of squares = 

641 9E-f-01 


Degrees of freedom = 10 



Par amet er 


es t 1 mat e 

Std error 

t “Value 

BC 0 ) 


114SEi-0a 

1 93£E+00 

.59£6E+0£ 

B( £) 


.S9a7E^00 

1932E+00 

. 3068E+01 

Bf 3) 


75£3£'*-00 

1 932E+00 

3894E+01 

BC 4) 

. 

4402E'f 00 

1932E+00 

££79E+01 

BC 6) 


749BE4-00 

1 93HE-1-00 

388£E'^01 

BC 7) 


70£3E ^00 

193£E-t^00 

.3636E-I-01 

BC 8) 


88£7£+ 00 

193£E^“00 

.4569E-i'01 

BC 10) 


.3E0£Et00 

1 93cE-t-00 

181 3E+01 

BC 11) 

- 

IITOE-^01 

1938E+00 

.6056E+01 





ANOVA 

FOR THE PRESENT 

MODEL 


,t i 

f f ^ t f t * 'f' t i* is -f It 5ft 

(: is i; fit 


Source 

DOF 

Sum Square 

Mean Square 

F“Va 1 ue 

MODEL 

9 

E49E+04 

. £76E+03 

. 3E0E+0 


FIT 

ERROR 


6. ,607E+01 .101E+01 .117E+0E 

4 346E + 00 .S64E--01 


1 9 . 


TOTAL 


. 249E+04 


. 1 31 E+03 


